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Abstract 


The concept of using an airborne laser doppler 
velocimeter is evaluated for diagnostics of the wind 
field associated with an isolated severe thunderstorm. 

Two scanning configurations have been identified, 
one a long-range (out to 10-20 km) roughly horizontal 
plane mode intended to allov; probing of the velocity 
field around the storm at the higher altitudes (4-10 km) . • 
The other is a shorter range (out to 1-3 km) mode in which 
a vertical or horizontal plane is scanned for velocity 
(and possibly turbulence) , and is intended for diagnostics 
of the lower altitude region below the storm- and in the 
out-flow region. 

The essential conclusion of the study is that air- 
craft flight velocities are high enough and severe storm 
lifetimes , are long enough that it is possible for a single 
airborne Doppler system, operating at a rar'^e of less than 
about 20 km, to view the storm area from two or more sub- 
stantially different aspects before the storm character- 
istics change appreciably. This permits resolution of 
the spatial distribution of the vt'locity vector by a 
single coaxial system operating in a fly-by mode. 

In addition, it is concluded that when either or 
both peak power rather than mean power and/or peak data 
handling rate rather than mean rate are limiting factors, 
the optimum signal-to-noise as well as optimum velocity 
reisolution are to be obtained by using pulse lengths com- 
parable to the overall dimensions of the viewed region 
and, as well, by operating in a continuous azimuthal scan 
mode. Range resolution is to be recovered in this mode 
by post-flight image reconstruction analysis. 
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CONCEPTUAL DESIGN OF AN AIRBORNE LASER DOPPLER VELOCIMETER 
SYSTEM FOR STUDYING WIND FIELDS ASSOCIATED WITH SEVERE 
LOCAL STORMS 

I. Introduction 

The purpose of this study is to develop design 
concepts for an airborne laser doppler velocimeter system 
for sensing the' vector wind fields associated with severe 
local storms. Presently available remote techniques are 
limited to microwave Doppler radar and give information 
primarily about regions of substantial radar reflectivity/ 
i.e., essentially/ v;here there is precipitation. The 
prime potential for a laser system is to regions of v/eak 
echo or no echo return for conventional radars. In general, 
this comprises three spatial volumes: the lov7 level 

boundary layer below the cloud structure, the internal 
updraft region (commonly is a weak echo region ) and the 
external velocity field surrounding the storm cloud, parti- 
cularly in conditions of strong vertical wind shear and/or 
for rotating storms. Convergence rates and shear or 
vorticity are of particular interest in the regions external 
to the cloud, v;hereas velocity and turbulence levels or 
mixing rates are relatively more important in the internal 
regions . 

A complete sequence of measurements must be made 
within the overall lifetime of the storm, typically 30 
minutes, v;ith sufficient time resolution to follow changes 
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in the structure of the storm (A,t ~ 2 to 5 minutes) . The 
coverage area and spatial resolution vary somewhat depending 
on which part of the storm is being diagnosed. The main 
body of the storm can be adequately diagnosed with a 
square coverage area of 10 to 20 km on a side. Spatial 
resolution of 200 to 500 meters is necessary to resolve 
the basic storm structure. Better resolution {£ 100 m) 
may be required for the more sophisticated measurements 
where convergence rates and vorticity are the desired 
quantities. A velocity resolution of the order of 1 m/sec 
is adequate to specify the internal wind fields and the 
low altitude gust regions. In the region external to the 
storm poorer spatial resolution may be permissible (per- 
haps up to 1 km) but better velocity and velocity gradient 
sensitivity is required to be able to measure' the conver- 
gence and vorticity fields. High spatial resolution is 
also required to discriminate regions of turbulence from 
non-turbulent but strongly sheared regions. 

This report is divided into three sections. Section 
ir (and Appendix II) describe the basic characteristics 
of severe storms and establish design criteria in terms 
of space-time-velocity resolution and coverage require- 
ments. In Section III various scanning configurations are 
discussed and tv;o of these are selected as candidates for 
further discussion in Sections IV and V. 
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The unique capability of the aircraft to view the 
storm area from several diff^,5fent aspect angles before 
substantial changes in the storm configuration take place 
allow deduction of the velocity vector as V7ell as permitting 
substantial improvement in the spatial resolution-velocity 
resolution tradeoff by post-flight processing. A detailed 
analysis and simulation of this effect is undertaken in 
Section IV and design tradeoff factors established relating 
the manner in. which instrumental range resolution and 
number of views can be exchanged to keep the same final 
spatial resolution- Finally, in Section V, the results 
are summarized and specific recommendations made for 
instrument design. 



II, Dimensional Characteristics of Severe Local Storms - 
Summary 

The basic morphology of a typical severe loc?il 
thunderstorm is shown in Figure 1. Features of particular 
interest to the designer of a measurement system include 
the magnitude and scale of the mean and fluctuating wind 
fields, the time scale for their change and the overall 
extent of the phenomena. In Appendix II v?e have attempted 
to summarize and classify these basic features for a variety 
of severe storm types. Included are single and super cell 
thunderstorms, tornadoes and funnel clouds as well as a 
variety of peripherally related phenomena; dust devils, 
gravity currents, etc. 

We can identify three main geometrically distinct 
regions associated with a typical storm. First, the main 
body of the storm, extending from 2 or 3 km altitude up 
to near the tropopause, is the primary visible object 
(clouds) and contains the regions strong up and dovmdrafts. 
It is generally opaque to optical sensors and laser velocity 
measurements V7ill be limited to the near surface regions. 

The second region is close to, the ground (up to 3 km) for 
a mature thunderstorm and is a horizontally extended region 
of cold air beneath (and ahead) of the main body of the 
storm. This region results from cold air in the storm 
sinking and spreading horizontally. Especially near the 
leading edge of- this outflow region, strong gusts may be 
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Figure 1(a): Typical Thunderstorm Cell Moving from Left 



Figure 1 (b) : Typical Thunderstorm Cell Moving from 

left to Right (top vie\v near ground) . 



encountered. The lov/er portion' of this region is typically 
clear air. Water droplet clouds or dust clouds may be 
present on occasion and occupy part of this region'. 

The third region characteristic of the storm is 
the inflow or outflow region surrounding the body of the 
storm at altitudes above the low altitude outflow region. 
The horizontal flow convergence and vertical shear of the 
wind field in this region describe characteristics of the 
energy and vorticity input to storm distinct from the 
vertical convection inputs. These may be particularly 
important in determining the tornado-producing potential 
of the storm. 

In Figure 2 we have summarized the essential 
•velocity-size characteristics of severe storms. 

Four essentially different phenomenology regions are 
identified. First,, the storm cell region itself comprises 
the main body of the updraft/downdraft region and is the 
primary visible feature of the storm. Because it is 
generally impenetrable optically, laser diagnostics 
are confined to measurements of the near surface velocity 
fields. 

Surrounding the storm is the external flow region. 
For single cell storms this is a generally clear air 
region of much 'lower velocity shear than the interior. 
Velocity measurements in this region look for the general 
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fluid dynamic structuire of the storm environment: the 

vertical shear, the vertical vorticity, and the convergence 
or divergence of the horizontal velocity. Three altitude 
intervals are of interest: 0 to 2 km (the cold out- 

flow region) , 2 to 8 km (the main vertical extent of the 
storm) and 8-12 km, the region of high altitude outflov? 
in the mature stage. 

Coverage dimensions are of the order of twice the 
storm crossection dimension in the 2 - 8 km interval 
(i.e. 5 to 10 km) up to several times these values in the 
low (and high) altitude intervals (i.e. from 10 to 50 km). 

The low altitude boundary layer region beneath 
and around the storm we define as the third distinct 
measurement region. This is the region of strong cold 
outflov; and may extend horizontally many times the storm 
cell dimension (10 - 50 km) although the vertical height 
remains limited (< 1 tp 3 km) . 

The fourth phenomenology region includes the small 
scale and localized features such as tornadoes, water 
spouts, and funnel clouds. 

In Figure 2 v/e have also indicated suggested spatial 
and velocity resolution design goals for these regions. 
These were established essentially as 10% of the charac- 
teristic dimension and characteristic velocity respectively 
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For the purpose- of defining a laser measurement 
system we have grouped the various phenomena into the 
first three general categories just discussed (excluding 
tornadoes) and their basic characteristics are presented 
in Table 1. As indicated in this table we have used these 
characteristics to define tv;o distinct measurement regimes, 
one appropriate to the low altitude outflov/ region and 
the other to the higher altitude region exterior to the 
central body of the storm. These are essentially the 
second and third regimes discussed above. The central 
body of the storm (the first regime) is judged to be too 
opaque for laser diagnostics except for cloud surface 
■measurements. These lat'ter v/e include as par-t of the 
exterior region. Measurement system candidates are dis- 
cussed in Sec-fcion III. 
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Characteristic 

Features 


Overall Horizontal 
Scale 


Vertical Scale 


Spatial Resolution 
Required 


Time Scale, . . , 

(lifetime) 


Time Resolution 
Required 



Low Altitude 
Outflow Region 


3-6. km 


- 10 km 


.100 - 300 m 


10 - 30 min 


6 .- 30 km 


.1,00. - 3 00 m 


.10 - 30 rain 


2-4 min 


2 - 4 min 


External 

Region 




- 10 km 


- 10 km 


■■;300 - 600 m 


10 - 30 rain 


- 8 min 


Typical Horizontal 
Velocity 

■ 

10 -■ ,2 0 m/s 

. 5 40 m/s 

' 

0.5-4 m/s 

Typical Vertical 
Velocity 

5-10 .m/s 

up to 10 m/s 

0.1 - 1 m/s 

Useful Measured 
Quantity 

6u^^ , w , u,.y , . 

. .... .U/ V,. 6u^ 

Vh , ?jxu , u,v 

Opacity to Optical 
System 

high, opaque 
except for surface 

Partially clear, range 
limited by atmospheric 
■ • absorption 

Clear for single cell 
storms, range limited 
by signal strength. 

Xaser Measurement 
Concept 

No 

Yes 

Yes 


Table 1. Characteristic Parameters for Three Storm Regions 




























Ill • Scanning Laser System Concepts 
III.l Introduction 

Two quantities are of basic interest; the wind 
velocity vector u. and, in regions of strong gradient, 
the turbulence levels (and', associated vzith these, the 
mixing rates) . In general the v/ind field has three vector 
components which vary from point to point. A single 
laser radar system operating in Doppler mode senses only 
one component and cannot by itself provide a direct 
general measurement of the instantaneous vector velocity 
of a wind field which is varying in an essentially random 
manner in space and time. Fortunately there are a number 
of good approximations that are often valid that will 
allov; meaningful measurements to be obtained even for a 
measurement system that cannot measure the general vector 
field. We group these into three categories: first, those 

that happen to be valid for the normal physical state of 
the atmosphere, second, those that depend upon some spatial 
uniformity of the wind field and, third, those that 
depend on a temporal uniformity. 

The dual Doppler diagnostic technique is an example 
of the first categor'y. Here two remotely located Doppler 
radars measure simultaneously tv/o of the three velocity 
components as a function of position throughout the 
measurement volume. The fact that v/ind velocities are 
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always small compared to the sound speed implies that the 
atmosphere is essentially incompressible and this condi- 
tion can be used to generate the third velocity component 
by integration of the continuity equation. The technique 
does require that the two-component velocity field be 
measured throughout the spatial volume and that a suitable 
boundary condition can be established (typically that the 
vertical- velocity vanishes at the ground surface). 

In the second category we include the well-known 
VAD technique (velocity-azimuth display) for measuring 
vertical wind shear. Here the vertical velocity as well 
as the horizontal gradients in the wind vector are assumed 
negligibly small so that measurements of the Doppler 
component at the same altitude but at different horizontal 
locations can be combined to yield the horizontal v;ind 
vector. 

Temporal homogeneity (third category) is used in 
the 'drift' type correlation measurements to deduce velocity 
vector components not directly detected by a Doppler system. 
Here inhomogeneities drifting across the line of sight 
are detected in some fashion (time correlation between 
signals from adjacent parallel lines of sight, or time 
spectral analysis of the response of a broad beam sensor 
whose aperture can be controlled so as to be selective 
for certain spatial scales.) , . 
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Almost by definition it is not easy to apply these 
techniques to a localized storm. The strong three 
dimensional spatial gradients of the overall storm V7ind 
field essentially preclude VAD concepts except perhaps 
on rather short and uninteresting length scales. Drift 
measurements require statistical homogeneity on a time 
scale many times the transit time associated \>7ith the 
correlation length scale. This type of technique is 
potentially applicable to a ground based storm monitoring 
system in which the optical sensor utilizes its very high 
azimuthal resolution to sense small scale inhomogeneities 
at each point in space. It requires a very high data 
aquisition load in its simplest version since the • spatial 
scale .sensed must be much smaller than 'the basic resolu- 
tion required. This requirement can possibly be mitigated 
by appropriate spatial filtering techniques which allov? • 
the density of small scale irregularities to be sensed 
without having to record their 'structure point by point. Even 
in this latter case, however, it is necessary to record and 
.analyze the data on a time scale very much shorter than 
the basic time scale associated with the oveirall V7ind field. 

It appears probable that this type of measurement 
will be most useful in an environment where properties 
are uniform in at least one direction so that the data 
recording rate is reduced to a reasonable level. 
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Airborne systems have difficulty in utilizing any 
type of measurement technique, such as this one, that 
require extended periods of time for each velocity 
measurement. The dwell time Or time period allowed for the 
total volxime measurement varies in proportion to the range 
in a fly-by mode and can only be made long at the expense 
of reduced signal strength and loss of spatial resolution 
associated with long range operation. 

The prime advantage of aircraft operation (in addition 
to the crucial one of mobility) is that, since the obser- 
vation platform moves at a speed very much higher than the 
mean speed of the storm, a single aircraft can obtain near 
simultaneous views of the storm area from several different 
aspect angles. Thus it potentially can simulate the capa- 
bility of several remote fixed location platforms. As will 
be shown, this capability can be converted into a capability 
for high spatial resolution in three dimensions. 

III. 2 Design Criteria 

We consider that the following problem areas, must be 
addressed in the design analysis: 

1. Ability to measure a) velocity vector 

and/or b) turbulence level 

2. Spatial Resolution 

3. Temporal Resolution and Dwell Times 

4. Propagation losses and signal-to-noise 
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5. Data Handling and Processing 

Each of these areas v/ill be discussed in the following 
sections. 


III. 2.1. Measu'reinent of the Vector Wind Field 'and Turbulenc e 
Field 

A Doppler measurement at a single space-time point 

can be made in a time a fev7 multiples of the reciprocal of 

5 0 

the Doppler shift (At ~ X/2V ~ ^(H/se ' c T whereas a 

drift or correlation measurement takes a time a few times 
L/v where L is the inhomogeneity scale. The minimum inhomo- 
geneity scale resolvable is determined by diffraction 
L ~ XR/D. At a range of 1 km with a 0.3 meter aperture, 

L - 0*03 meters and At "" 0 . 03y^v (m/sec) seconds. Thus for 
an area measurement v/ith 30 x 30 grid the net time required 
(for V = 1 m/sec) is 27 seconds for the drift system as 
compared to 0.004 seconds v/ith the Doppler system. Since 
times of the order of 100 times these values are required 
for a single measurement (to get 10% accuracy), Doppler 
techniques appear to be the only method applicable for a 
system v;hich m.ust obtain volume data in two or more dimen- 
sions in a serial processing -raode^ 
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In a parallel processing mode (each space point 
monitored continuously and simultaneously) drift type measiire- 
ments are conceptually applicable. Although v/e feel that’ 
it may be possible to develop such a concept in the future, 
no such scheme has presently been applied in a scanning mode. 
In any case, because of the small scale of the inhomogeneities 
sensed, the scene registration problems for an airborne system 
would be formidable in addition to the fact that observation 
times are necessarily short. For the present study it v^as 
felt that such a concept would require^ a substantial effort 
even to be defined, and detailed evaluation could not be 
undertaken within the present study framework. Thus it 
was concluded that only a Doppler' technique offers promise 
for early implementation. 

In the pure Doppler mode it is mandatory to have at 
least two different aspect angles to obtain the velocity 
vector. At a range of 10 km and a flight velocity of 150 
meters per second the viewing aspect angle changes at the 
rate of about one radian per minute. Since the major char- 
acteristics of the wind field change over substantially 
longer time periods (see Table 1) , a single airborne moving 
platform can, in principle, be used to measure the velocity 
vector. The error in one component due to the- fact that 
the components are measured at different times is of order 
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s ■ R ' 

AV ^ ^ v/here is a typical storm velocity^ the 
s 

aircraft speed, R the range and the characteristic time 
associated with the development of the storm wind field. 
For R = 10 km, T = 10 minutes, and V_ = 150 m/sec , 



Of course, this error can be substantially further reduced 
by averaging one set of measurements v/ith a second set taken 
during an immediately follov;ing fly-by on the same course but 
in the reverse direction. Here the errors to first order 
have opposite signs. 

The turbulence field can be measured in tv7o v;ays: 
a) from the width of the Doppler spectrum and b) from the 
fine scale structure of the measured vector velocity field. 

In principle both can be measured as scalar quantities from 
only one aspect angle- However, the V7idth measurement V7ill 
•be contaminated with spectral broadening effects associated 
with v/ind shear. The spectral broadening (in velocity units) 
due to wind shear along the line of sight is 




AV 


U 


9S 


r AS 


v/here AS is the longitudinal resolution. In a simple shear 
layer, the rms v;idth due to turbulence is 
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AV 


rms 





where Z is the largest eddy scale length . (mixing length) 
contributing to the random turbulent field and is the 

velocity shear driving the turbulence. Mixing lengths of the 
order of 10 or 20% of the thickness of the shear layer 
are expected. Since the requirements of the measurement 
system necessitate a spatial resolution AS also of this 
order or smaller (i.e. AS < 0.1 L :: z) the spectral broadening 
due to turbulence and that due to longtitudinal gradients 
in the mean v;ind field may be comparable. Thus we- conclude 
that the spectral V7idth of the Doppler return is not likely 
to be an unambiguous estimator of turbulence level and it 
will be necessary to measure directly the local fine scale 
irregular structure of the velocity field. Hovzever ^’•e 

the small scale turbulence is expected to be isotr , iC/ xt 
would be possible to obtain an adequate measurement of 
the turbulence level from a single viewing aspect. The 
requirements on range resolution for such a measurement are 
higher than for the mean velocity field- by perhaps a factor 
of 3 to 6 Ci-e^ ^^here a range resolution of 300 meters may 
be adequate for the mean velocity fields 50 to 100 meters 
vzould be required for a turbulence measurement. 
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11,1.2.2. Spatial Resolution • 

The required spatial resolutions for various regions * ; 
are defined in Table 1. Essentially a resolution at least 
as good as 100 meters is specified for the closer ranges and, 

V 

at all ranges, for turbulence measurements. At longer 
ranges the requirements are weaker for the mean velocity 
field (300 meters) . 

The CW laser systems obtain range resolution by focus- 
ing, and the effective length of the resolved volume increases 
strongly with range (see Figure 3). Reference to this 
figure shov7S that, at least v;ithout additional processing, 

CW systems fail to achieve adequate resolution at ranges 
■greater than about a kilometer (assuming a useful effective 
aperture of 0.3 meter diameter). The actual useful aperture 
achievable depends, in addition to any mechanical or cost 
constraints, on both the range and the ambient turbulence 
level along the line of sight. Limiting effective apertures 
can be determined from Figure 4 . 

Pulsed systems achieve range resolution directly 
by the finite pulse duration. Associated with the finite 
pulse length is a finite frequency bandwidth and velocity 
uncertainty.- In Figure 5 we shov; how the jivailable range 
resolution- velocity resolution product compares to the 
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typical turbialence levels (see also Fig. 2 ) . It is appar- 
ent that coarse velocity information about the larger scale 
features can be obtained directly v;ith a pulsed system V7ith- 
out additional processing. • Improved range resolution will 
be needed however for accurate measurements. 

There are tv;o methods that potentially can be 
implemented to improve this resolution: range processing 

and aspect processing. 

Range Processing 

When the signal-to-noise ratio is high, data taken 
as a function' of range can be processed so as to improve the 

effective range resolution. The improvement achievable 
depends on both shape of the range response function and the 
signal-to-noise ratio. The optimal method for accomplishing 
the improvement in the presence of white noise is called 
matched filtering and is described in detail in .Appendix L . . 
Substantial improvement at high values of SNR can be 
achieved only V7ith range response functions that have sharp 
cutoffs. Fuzzy or spreadout functions permit only modest 
improvements. 

In Appendix I it is shown that the optimally processed 
range resolution AX* may be approximately estimated at high 
SNR from the expression 

|g(k = 2/AX*)/g(k> 0) ] = 


23 


where SNR is the ( power ) signal-to-noise level and g(k) 
is the spatial Fourier transform of the range response 
function. For a CW focused system the range response func- 
tion at ranges where good depth of field can be achieved is 


^ . constant • 

f (x) = , - 5 — 

Here AX is the range resolution betvjeen 3 db points. In 

CW 

Fourier space this resolution function is 

gCh) ?= constant » exp j k.J AX^^/21 
Thus, the ‘best* processed range resolution is estimated as 

AX / 

AX* = AX „ An (SNR An (SNR ) 

cw/ voltage 2 / pov;er 

Power signal-to-nbise ratios betv/een 20 and 30 db permit 
range resolution improvements by a factor between 2 and 4. 

The same type of argument applied to the improvement 
of the range resolution of pulsed systems. For a perfectly 
square pulse 

f(X)-= 1 for |x| £ -“^pulse ’ 

= • 0 otherv/ise 
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The Fourier transform is 


sin ((TT/4) k. ‘Ax_^^ 

9 0 ^) = — 


se 


The optimal range resolution is obtained from the relation 


• /sin(l. 




SNR 


,pov;er 


x.e. 


AX* AX T ^/ 2 /SNR 

pulse V pov;er 


(The factor of 2 came from averaging the sin squared term. 


This, expression indicates that substantial improve- 
ment is possible for such pulses at high signal-* to -noise 
ratios Cfactors of 10 to 20 for values of 20 to 30 db) 

However, the result at high SNR is totally dependent on the 
sharpness of the leading ; and trailing edges of the pulse. 

Finite rise and fall times of the actual pulse will strongly 
suppress this dependence - 

In summary it appears that, at high signal-to-noise 
ratios, significant improvement in range resolution can be 
achieved by processing the data in range for both CW 
and pulsed systems. The improvement anticipated at an 
SNR = 30 db varies betv;een a factor of about 3 to an upper 
limit of about. 20 , The higher values require very precise 


shaping of the range resolution function and are undoubtedly 
difficult and expensive to achieve in practice. Nevertheless 
some improvement can always be achieved if the SNR is high 
enough and such processing should be included in the design. 
The optimal (as v/ell as the most- direct) technique is that 
of matched filtering. 

Aspect Processing 

All optical systems basically exhibit very high 
lateral or' angular resolution but much poorer range resolu- 
tion. However, by making observations at various aspects 
the range resolution may be vastly improved. A hologram 
is an interesting example of this process. Whereas a . 
normal photograph of a scene essentially records only the 
lateral variation and has virtually no range resolving 
poxv^ef , a hologram ^is able to record the scene as it appears 
at a variety of aspect angles. This record can be viewed 
subsequently using various methods to extract range infor- 
mation {triangulation, parallax, depth of focus). 

. A technique which might be described as a type of 
incoherent holography is currently finding extensive appli- 
cation in a variety of fields, particularly in medicine in 
v/hole body X-ray scanning. Basically a number of one- 
dimensional projections of a two-dimensional, (crosssectional) 
transparent object are formed at a number of equally spaced 
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aspect angles. These projections are then combined together 
by a computer to yield a picture of the two-dimensional 
structure. In medical applications the number of vievzs 
required for useful reconstruction ranges from 2Q to 200. 

An example of such reconstruction applied to x-ray imaging 
of the cross-section of a human body is shov/n in Figure 6- 
A similar technique is potentially applicable to the recon- 
struction of velocity fields when the scanning measurements 
can be made at a number of yiev/ing aspect angles . A detailed 
analysis of the mathematics of this procedure is presented 
in Section IV. For the discussion here, hov/ever, it is 
sufficient to note that significant image improvement is 
possible with as few as four views although a much larger 
number may be desireable to recover the fine scale detail. 

An attempt v;as made in the analysis of this method in 
Section IV to establish a quantitative relation between the 
spatial resolution improvement and the number of aspect 
views. Because a matched filter analysis was used the result 
depends explicitly on the signal-to-noise ratio. The computer 
simulation technique used and discussed in Section IV yielded 
satisfactory results for a scalar field (such as the turbu- 
lence level) and roughly exhibited an inverse relationship 
between number of vievzs and processed spatial resolution 
(i.e. AX* ” AX‘^/N^, where A is the single sensor resolution 




which here contain* an X-ray opaque ilye for diajnortic purposes. 
The two roundish areas to the left and right of the spinal column 


CROSS SECTION OF ABDOMEN of a living human subject ap- 
pears in this X-ray picture made by reronstrvetion from projec- 


tions. This section is seen from the same perspective as the one on 
the cover of this issue: from above the subject’s head. The vertical 
bar at the right is the key to the X-ray density of the tissues in the 
section; white represents the greatest density, black the least den- 


are the kidneys. The small dark areas extending from the top 
around to the lower left are sections through loops of intestine. 
This picture, like the one on the cover, was made by the Delta 
Scanner built by Ohio-Nuclear in the course of a study by Ralph 


sity and the colors the densities in between. The white .area* below 
the center of the section are the spinal column. The large light- 
colored area running from upper right to lower right is the liver. 
The light area near the upper end of the liver is the gall bladder. 


J. Alfidi, M.D., of the Cleveland Clinic Foundation. The picture 
required 180 separate X-ray projections, made at one-degree inter- 
vals around a half<ircle. Total dose of X rays was approximately 
the same as that needed to make a single conventional X-ray plate. 


Figure 6. X-ray imaging of the cross-section of a human 

body. (Scientific American, Oct. 1975, Page 57) 
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and the mimber of views) . The corresponding dependence 
for the vector v/irid field v/as less simple, and it is felt 
that further analysis will be required to, establish an 
optimum processing technique. The tradeoff betv/een number 
of views and spatial resolution (derived from the present 
simulations) is shown in Figures 24 to 26. 


111. 2. 3 Temporal Resolution and Disrell Times 

The dwell time available to examine each resolu- 
tion’ volume may be estimated as follows. The time 
available to examine one resolution volume is the total 
observation time divided by the number of resolution 
elements. To scan a single plane this is 

At = ^ 

N N 
a X y 

where R is the mean range and the flight velocity. For the 
two scanning modes discussed in Section III. 3 (50 resolu- 
tion elements in range, 25 in azimuth) this dwell time is 
100 milliseconds (for a maximum range of 20 km) or 25 . 

milliseconds (maximum range of 5 km) . The maximum per- 
missible integration time is this time increment divided 
by the number of views. However, vrhen post -flight process- 
ing is carried out an incoherent integration is automatically 
performed (albeit V7ith some losses) . Thus, the signal-to- 
noise ratio can be estimated based on the above net dwell 
times if some account is taken of the processing losses 
(such as scene registration losses) . The maximum signal- 
to-noise is to be obtained v;ith the highest mean output power 
together with the lowest prf consistent V7ith the data 
taking mode. However, these values of the prf are generally 
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quite STDall^ For a 100 meter resolution and 15 0 m/s 
flight velocity^ the minimum time between pulses should 
be ■0.7/N^ sec. v;here is the number of viev7s. For = 

10 -the optimum prf is only 14/sec corresponding to a 
pulse ehergy of 300 millijoules at a. mean power of 5 watts. 
At lower pulse energies and higher prf *s some losses 
are to be expected, due to the fact that incoherent rather 
than coherent integration takes place. 
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111.2*4. Propagation Losses and Signal-to-Noise 

Atmospheric losses due to molecular absorption 
and turbulence are significant at lov; altitudes at the 
range of interest even under good seeing conditions, and 
limit the usable ranges. 

For a pulsed system operating in a collimated mode 
(focused at infinity) , the SNR is 


S/N = n 


D no- ct 


32 hc/A 


R' 


(^) 


L(R) 


where 


J = . net output energy in the observation time 
D = aperture diameter 
T = pulse time 

■na - particle density times scattering cross-section 
g =r (mean value at range R) 

L = detection efficiency 
R = range 

L = propagation loss factor 

The factor L is the product of the molecular loss 

and the turbulent loss L. , . 

turb 

At ranges less than the Rayleigh distance (R._,. = 

2 ■ . . 

TfD /4A) , the signal-to-noise ratio vanes with range only 
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as a result of the loss factor f«; 


s/N - n — ^ ^ ■ l(r). for R < R* . 

2tt hv D^. 

For a 30 cm effective aperture diameter and A = 10. 6y, 
this expression is appropriate for ranges less than about 

7 km. 

Within the range limits, some benefit can be 
obtained by focusing on the range interval of primary 

interest. For pulse lengths less than the focal depth 

( ct 2 2\ 

< 4R X/D ) , the signal-to-noise from the region near 
the focal point is 

2 

. 0"X b nc CT T^■0^ 

32K5 

and exhibits an inverse square dependence on range in 
addition to the range dependence through the propagation 
loss factor. 

Molecular absorption results principally from C02‘ 
lines and the wings of the nearby water vapor bands . 
Signal-to-noise estimates for a range resolution equal to 
10% of the actual range are given in Figure 7. The 
signal-to-noise ratio is given in terms of its value per 
unit energy in a resolution volume. For coherent inte- 
gration of successive pulses the actual signal-to-noise 
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for in^ny pulses can he obtained directly by multiplying 
the. ordinate by the total laser energy allocated to a 
given resolution volume. Only pulses separated in time 
by less than the laser line correlation time (reciprocal 
of the laser bandwidth) can be integrated coherently, 
however. For the more usual case where the time between 
pulses is greater than the correlation time, incoherent 
integration must be used and the effective signal-to-noise 
will be reduced from the coherent value .given in Figure 7 
by the facto.r N v;here N is the number of incoherent 

pulses integrated. 

The sxgnal-to-noise levels given in Figure 7 are 
based on a back scattering coefficient of 10 km at 
sea level (which is approp.riate for clear air) and an' 
overall instrume.nt-detector efficiency of 0.6%. Actual 
SNR values may vary upwards substantially from this value 
in hazy conditions and/or for improved sensors. Neverthe- 
less the basic conclusion that must be drawn from this 
figure is that molecular absorption will limit the 
maximum vrarkable range to some value less than about- 5 km 
at sea level. 

Under conditions of st.T-\,ng surface heating, low 
altitude turbulence will furth. degrade the returns. 

The signal-to-noise ratio in this case is conveniently 
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signal/noise per joule (77= 1.0) 


L-’.icjuro 1 , " Signal-to-Noise per Joule vs. Range 

(At-.mospheric ■ attenuation included). 



'signal/noise per joule (19 = 0 r 006 ) 


expressed in terms of an effective coherent aperture 
radius r^. When r^ is small compared to the actual aper- 
ture D,/ the optimum signal-to-noise (that obtained v/hen 
. the beam is focused on the range point of interest) is 
given by 

S/N = aina^cT . 

K 

The value of r is given by the expression 

a 



where R is the range and the refractive index structure 
coefficient due to the turbulence. 


Ill, 2.5, Datci Handling and Processing 

The average rate at which data is to be acquired . 
and recorded is relatively lox^?. For example, for a planar 
scan with a 25 x 50 mesh, allowing 2 minutes for the total 
observation and 10 separate scans or viev7s of the observation 

plane the mean data taking/recording rate is about 200 

values/second. With a pulsed system, however, the peak 

data rates are very much higher; the Doppler spectrum 

must be analyzed within one pulse time (= 2AX/C) . For a 

range resolution of 300 meters this is 2.1 ysec (0.7 ysec 

for 100 m resolution) . This corresponds to a peak data 

5 5 

handling rate of betv;een 5 x 10 and 1.5 x 10 output 
values/sec. Tv;o techniques are possible to reduce, this 
peak load: either buffering the input or using longer pulses 

but maintaining the high spatial resolution by post-flight 
processing of the multiview data. In the design a tradeoff 
analysis should be performed to compare the advantages of 
lower peak data handling rates to the increased complexity 
of the multiview process. It should be also noted that 
longer pulses also improve the coherent to incoherent inte- 
gration time ratio. 

When the signal-to-noise ratio per pulse is signi- 
ficantly greater than unity, optimum use of the available 
power requires that the signal energy be spread as uniformly 


as possible over the individual measurement degrees of free- 
dom. For a tv70-dimensional field of the type we are con- 
sidering here , it is convenient to examine the V7ay infor- 
mation is collected in the (spatial) Fourier transform 
plane rather than physical plane. Here the coordinates 
are the v;avenumbers and and the measured data consists 
of the amplitudes of the Fourier modes v;ithin a circle 
about the origin having a radius k = ir/Ax, where- Ax is 
the limiting spatial resolution. The spacing betv7een 
different or independent modes is Ak = '2 tt/L where L is 
the overall linear dimension of the scene. The objective 
of the measurement is to sample and measure each of these 
modes . 


As is discussed in some detail in Section IV the 
effect of low range resolution in the laser spatial re- 
sponse function is to permit sampling of the k , k space 

X y 

by a single sensor only in the neighborhood of a line 
emanating from the origin (k^ = k^ = 0) and perpendicular 
to the direction of view- Full sampling requires super- 
position of several viev7s each at a different aspect angle 
(see Figure 11) . 

In the presence of white (instrumental) noise opti- 
mum measurement strategy requires that the signal energy 
should be spread more of less uniformly throughout the 

■y 

sampling area in k space. Emphasis on particular directions 
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at the expense of others is to be avoided. For a multiview 
system with a finite number of views (say 10 viev/s) it is, 
therefore, better to have 10 different aspect views each 
having a signal-to-noise ratio of 5 (=7 db) than only 5 
views each v/ith twice the signal strength S/N =10 (10 db) . 
This effect is shown for example, in the simulations of 
Figures 16 and 17. Here the same total signal energy was 
assumed constant but distributed over a varying number of 
view angles. The one-viev/ data in Figures 16 and 17 were 
calculated as the average of 8 identical measurements 
each with an independent noise realization; similarly the 
8-view data were obtained with a different noise realiza- 
tion for each view. It is clear from this figure that 
substantially more information is obtained v;hen the meas- 
urements are distributed uniformly in aspect angle. 

Thus the basic design strategy should be the following: 
After the desired spatial resolution (Ax) and coverage 
dimension (L) are established, the merits of v7orking with 
single sensors having poorer range resolution (than Ax) 
are to be assessed with respect to data rates, peak pov;ers, 
signal energy, velocity resolution, etc. If a substan- 
tialiy improved performance can be obtained in this manner, 
a multiviev;. concept should be generated. The minimum 


angular spacing betv;een independent viev/s is A0 = 0.5 


max 


Ax/L. 
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If this spacing cannot be achieved then there v/ill be a 
penalty in terms of the ultimate resolution achievable. 

The consequences of such incomplete sampling and attempts 
to ameliorate it are discussed at some length in Section 

IV. 

No attempt will be made here to discuss various 
specific data handling problems, such as storage, indexing, 
real time vs. post'-f light processing, scene registration, 
scene overlap losses. Required data recording rates for 
particular scanning concepts are discussed in the following 
section. 

Curves relating hoxv the number of separate views 
required must be increased as the pulse time is lengthened 
in order to maintain a given resolution are given in 
Figures 24, 25 and 26. These curve s.v7ere derived by analyzing 
the results of a.’ limited number of computer simulations 
of the multiview analysis. The results .depend on-.signal- 
to-noise levels as well as on the filtering processes 
assumed. The processing techniques used were not optimum 
and, for actual design purposes, it would be desirable to 
investigate improved raethocis, particularly for the vector 
v/ind field analyses. 
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III. 3 Selected Design Concepts 


III., 3.1 Introduction 

Based on the preceding analysis we have selected 
three separate design concepts for more detailed discussion 
Tv 70 of these/ a short range and a long range single plane 
scanning configurations are recommended for more detailed 
design. The third configuration requires viewing of 
several parallel planes simultaneously and consequently 
has a much higher mean data rate and more stringent require 
ments for signal-to-noise levels and scene registration 
requirements. It is suggested that detailed consideration 
of this last configuration be premised, on the successful 
demonstration of the simpler modes. 

III. 3. 2' A 'Short Range System 

This concept is intended to be appropriate for 
examining the low altitude outflov? region or for probing 
the vertical velocity profile v^ithin 1 to 3' km of the 
aircraft. In the simplest form the laser line of sight 
makes one of a set of fixed angles with respect to the 
aircraft flight axis. A coverage area whose dimension 
is smaller than the mean range is contemplated. As the 
aircraft completes a sv/eep of the beam through the 
coverage area^ the beam aspect angle is reset so as to 
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obtain an additional sweep. This procedure is repeated 
for as many different aspect angles as geometry permits. • 
The configurations for a tv70-aspect and a three-aspect 
angle, observation are indicated in Figure 7a. In prac- 
tice for coverage of an area v;hose lateral dimension is com 
parable to the range, only tv/o or three aspect angles are 
possible. In Table 2a we give characteristic values for 
the parameters for such a configuration. 

Better area coverage and more aspect angles can be 
obtained at the expense of a more complicated optical- 
mechanical system in a continuous scanning mode (see Figure 
7b) . Here a full scan in azimuth is made in the time it 
takes the aircraft to m.ove one resolution length. The 
mean data rate is much higher than in the fixed position 
line of sight (LOS) configurations. Suggested parameters 
for such a scanning configuration are given in Table 2b. 
Advantages of the continuous scan multiview mode versus 
the 2 or 3 position LOS mode, include better area coverage 
and smaller scene overlap losses, lower peak power require- 
ments, lower peak data handling levels , lower incoherent 
integration losses, good velocity resolution. Disadvantages 
are the more complex scanning requirements, the requirements 
for additional post-processing of the data, and the (poss- 
ible) requirements for greater data storage. If no in- 
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Table 3a 


Suggested Characteristics for a Two-and Three-View 
Fixed Position LOS Side Looking System. 


Assumed Flight Velocity (m/s) 

Range to Center of Coverage Area (km) 
Diameter of Coverage Area 
Spatial Resolution (meters) 

Useful Grid 

Total Observation Time 

Kean Time Between Azimuthal Points 

Mea.n Data Recording Rate (value/sec) 

Number of Data Points per Frame 

Viev; Angles 

Kean Ranges in Each View (km) 

Observation Time in Each View 

Peak Input Data Rate (value/sec) 

Pulse Length (time) 

Pulse Energy/average watt of 
Laser Power 

Optimum prf (sec ^) 

Peak Pulse Power at Optimum 
prf/v/att of Laser Pov/er 


Short Range System 
2 Viev7 3 View 

150 m/s 

3 km' 

3 km 

125x125 m 
25x25 
47 sec 


0 . 9 sec 
-1 

28 sec 
1320 

60 °, 120 ° 


2 . 5 km 
125x125 m 
20x20 
63 sec 
1.0 sec 
20 sec ^ 
1260 

45°, 90°, 135° 


Incoherent Integration loss at prf=100sec "" -10 db 

at prf=1000sec ^ -15 db 


3 . 5 , 3 . 5 km 4 . 2 , 3 : 0 , 4 . 2 km 
23,23 sec 23,17,23 sec 
1.15x10^ sec''^ 

125 m (0.4 ysec) 

0.9j/watt 

1.1 sec ^ 

2000 KW/w'att 

-1 


Long Range System 
2 View 3 Viev7 

300 m/s 

15 km 

15 km 

300x300 m 

50x50 

2 min 


1.2 sec 
-1 

42 sec 
5040 

60°, 120° 


12.5 km 
300x300m 
42x42 
2.7 min 
1.3 sec 
32 sec ^ 
5180 

45°, 90°, 135° 


17.5,17.5 km 21,15,21 km 
1 . 0 , 1 . 0 min 1 . 0 , 0 . 7 , 1 . 0 min 
5x10^ sec ^ 

300 m (1 ysec) 

1. 2 j /watt 

0.8 sec ^ 

1200 KW/watt 

-10 db 
-15 db 


Table 3b. Suggested Charactr-ristics for a Multiview Continuous 
Scanning Side Looking System 


Assumed Flight Velocity (m/s) 

Mean Range (km) 

Spatial Resolution (meter) 

Frame Dimensions 
Useful Grid 
Number of Z^zimuth Points per Scan •' 

Useful Number of Range Points (?s total number) 
Number of Views per Frame 
Time Between Scans (sec) 

Mean Data Recording Rate (values/sec) 
Number of Data Points per frame 
Total Time for One Frame 
Peak Input Data Rate 
Pulse Length (time) 


Short Range System 
150 m/s ' 

3 km 
125x125 m 
3x3 km 
25x25 
25 
25 
30 


.85 sec 

700 sec ^ 

18750 

22 sec 

.4 


Pulse Energy /Average Watt of Laser 
Pov7er (joules) • 

Energy Per Processed Data Point 
Optimum prf (sec ^) 

Peak Pulse Power at Optimum prf/watt 
of Laser Pov;er 


Incoherent Integration loss at a prf=100 sec 

=1000 sec 


5x10 
3 km (10 ysec) 
0.03j 

0.9j 

30 sec"^ 

3 KW/watt 

-1 


Long Range System 
300 m/s 
15 km 
300x300 m 
15x15 km 
50x50 
50 
50 
50 

1.0 sec 
2500 sec ^ 
125000 

50 sec 
10 ^ 

15 km (50 ysec) 
0.02j 


l.Oj 


.-1 


-1 


-2.6 db. 
-7.6 db 


50 sec 
0.4 KW/watt 

-1.5 db 
-6.5 db 


flight processing is carried out, the data recording rate 
for the continuous scanning multivievr mode is 10 to 20 times 
greater than for the simpler tv;o view fixed LOS mode and can 
lead to a substantial data load. Hov/ever, the post-process- 
ing techniques are relatively straightforward and it 
should be possible to carry out in-situ and perhaps even 
real-time processing. Basically 1 or 2 two-dimensional 
Fast Fourier transforms are required for each azimuthal 
scan {i.e. roughly each second) with a 32 x 32 or 64 x 64 
grid. When this can be done the signal data recording 
requirements are relatively low (less than 100 values/ 
second) . 

111. 3. 3 A Long Range System 

The simplest version of the long range system has 
the same parameter structure as the short range system. 
Here, however, the range resolution requirements are less 
stringent. Candidate sets of system parameters are given 
in Table 2. Except for the differences in values of the 
system parameters, the configuration s are the same as for 
the short range systems. 

111. 3. 4 Multiple P l ane Scanning System 

In Figure 8 v;e shov; the basic configuration for a 
long range scanning system. This system would have the 
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capability of scanning and recording the Doppler signal 
from set of horizontal planes stacked in the vertical 
direction. Multiple looks at different aspect angles 
yield the components of the velocity in each scan plane. 

The already developed techniques for processing 
this two-component data to yield the third velocity 
component can be implemented. Here the vertical velocity 
is calculated from' the incompressibility relation 

^ ^ _ ly 

ax 3y 

The vertical velocity w may then be evaluated from 



In practice a more complicated equation must be solved 
to account for the fact that the individual scanned 
planes are not all parallel (see Fig. 8) . Here the 
measured in-plane components are related to the actual 
horizontal (u,v) and vertical (w) velocities according to 

Su =u'cos0 + V7 sino 

= VCOS0 -f- w sino 

The vertical velocity is then to be evaliiated from 


/ 


8w 

dz 


= tan 9 


/ 9v; 


9yJ 


^ /3s as 

1 r -h 

9y 


COS0 \, 9x 


No detailed analysis of this mode has been under- 
taken during this study. The mean data rates are an order 
of magnitude greater than for the single plana system 
and the requirements for measurement accuracy are greater 

Cin order to measure the gradients) . It v^as felt that a more 
definitive analysis of the simpler modes should precede 
analysis of this much more sophisticated technique (see 
Table 3 for tv7o possible configurations) . 
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10- 20 km 


Figure 8 Scanning Config 

a Long Range Multipi 
Scanning System 




Table 5* Possible Parameters for a Multiple Plane Scanning 


System 

Assumed aircraft speed 300 m/sec 100 m/s 

Mean range to coverage area 18 km 18 km 

Diameter of coverage area 15 km 15 kra 

Total observation time (2 viev;s ) 2 minutes 6 min. 

(3 views) 3 minutes 9 min. 

Spatial resolution 300 m x 300 m 

Mean data input rate for 10 planes 420/sec 140 sec 


Time between Vertical Scans 


1 sec 


3 


sec. 


IV. MULTI VIEW SYSTEM ANALYSIS 
IV. 1 Introduction 

We treat the LDV scanning system as an imaging sensor 
which views a scene from various directions. The scene is 
characterized by a reflected or scattered light intensity 
distribution a (x) and its Doppler spectrum. The imaging 
system records a degraded representation of this scene 
(an N-dimensional scene in general is recorded as an N- 
dimensional pattern v/ith finite spatial resolution) . The 
finite resolution or spatial gain of the recording system - 
is described by a response function G(x^-x,r). In this 
section we v^ill develop a mathematical analysis of this 
type of sensing system and the processing required to 
reconstruct the original scene. 

The recording system depends on two parameters: 
the viewing direction defined by the unit vector r and the 
displacement (x^-x) of a spatial point x' from the point 
X on v/hich the instrument . is focused. The geometry is 
shown schematically in Fig. 9. The 'photographic' mode 
is obtained when the resolution in the direction parallel 

A 

to r is so poor as to result effectively in a tv/o rather 
than three (or a one rather than tv7o) dimensional- record. 

In the presen-t analysis we will be concerned only v/ith 
sensing systems that approximate the mode for v/hich the 
resolution in range is very much poorer than the resolu-hion 



in azimuth. Good spatial resolution is then achieved only 
perpendicular to the viewing direction. Depth perception 
or range resolution is achieved by invoking the stereo 
effect (i.e. triangulation or parallax) . To utilize this 
effect to obtain depth resolution, comparisons and correla- 
tions must be made bet\r?een the various projections obtained 
by viewing in more than one direction. In this V7ay the 
good lateral resolution is used also to achieve localization 
in range (depth) . In the normal photographic viewing sit- 
uation the target consists of discrete points which can be 
identified and correlated in the different viev?s. However, 
in a diffuse scene, there are no such sharply defined 
features and the ability of an observer to sense depth from 
even two widely separated views deteriorates. 

I 

It is the purpose of this analysis to construct a 
model of the multi-viewing process for generalized scenes and 
to construct a mathematical procedure for combining the 
various projections into a processed full dimensional scene. 

Normally a three-dimensional scene is viewed and 
recorded as tv7o (or more) two-dimensional projections 
(pictures) . For the present analysis however we restrict 
ourselves to sti'ictly two-dimensional scenes v?hich are 
recorded as a number of imperfectly resolved two-dimensional 
scenes or, in the ’photographic* mode, as one -dimensional 
projections. 
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We will first (Sec. IV. 2.1) develop the theory 
for the imaging of a scalar field, such as a variable 
density of scatterers cr(:x,y). Secondly v/e will consider 
in Sec, IV. 2. 2 the. imaging of a vector field i.e. typically 
a wind velocity vector field by extending the analysis 
developed in Sec. IV. 2.1, 

IV. 2 Theory ' 

IV. 2.1 SCALAR MODE 

We consider the record of a two-dimensional scene 
as observed vrith a recording ins triment viewing the scene 

• A 

in the direction r with finite spatial resolution. When 

the instrument is set to measure the scene brightness 
">• 

at the point x, the actiial contribution to the recorded 
value at this point comes from points in the scene in 
the neighborhood of x: 

CO CO 

S(x,r) = f f G(s,r) a(xts) ds ds (1) 

j j X Y 

wt- CO CO 

•y ^ . 

Here G(s,r) is the resolution function of the instrument 
and cr(x) is the local brightness of the scene. The unit 

A 

vector r represents the viewing direction. We are 

-y 

particularly interested in resolution function G(s,r) 
of the general forra shox'^n in Fig. 10 in v/hich the lateral 
resolution is very rauch better than the longitudinal • 
resolution. 
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Fig. 10. Contours for the resolution function G{s,r) 

(see Eg. 5 for definition of and A^) 

The spatial Fourier Transform of Eg. 1 has the 

form 

s(k,r) = g(k^r) a (k) (2) 

Conceptually it is possible to recover the scene properties 
(i.e., .'cr(x) or, eguivalently , cr'(k)) by dividing by g(k,r): 

a'(S) = 2i|4l (3) , 

g(k,r)- 

However, in the presence of finite noise or measurement 
error, this result is very poorly evaluated in regions 
V7here g(k,r) is small. Here the noise (n) contained in 
the measured signal • 

s (k,r) = (1<,r) + n ,r) (4) 
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is strongly amplified and could destroy the image. 

The objective of the stereo viewer now is to 
utilize the fact that data has been taken at different 

/s 

values of r. Thus he wants to choose r as a function 
of k(=r^(k)) such that g(k,r^(k)) is never small. In 
this way different Fourier components are to be evaluated 
from different views. 

For example, consider the resolution function 

(5) 



This is an asymiiietric Gaussion v/ith its long axis (taking 
pointed in the direction r. The transform 
g(k,r) can be explicitly evaluated: 


g(k,r) = 


OO CO 

J I 

CO - CO 


G(s,r) ds ds 

X y 


f ->''2 2 ->^2 2 

= exp i - (k-r) A, /4 (k^r) /4 


} 


(6) 


For good recovery of the object scene v;hen A^ << A^^ 
the solution cr(k) should be evaluated from Eg. 3 with 
k*r maintained as small as possible, i.e., if possible 

near k “ h k V7here n is a unit vector perpendicular to 




simulated Example 


In order to demonstrate the viewing process, we 
consider a simple distribution of diffuse targets 
described by: 


a (X) = ^ “n A’ ) 


and v/e imagine that the observer views this field with 
the response function: 


G (s,r) 




(s*r) ^ ' 


(sxr)^ 


For values of a much smaller than At and and 

n 1 2 

as v;ell less than the mutual separations, the. targets 
appear as individual discrete points. On the other hand, 
when the are large, the scene brightness is smooth 

and slowly varying. 

The measured signal in the absence of noise has 

the form 


S^(x,r) 


I 


^CO —CO 


2 

(s- r) 


(sxr) 


(x+s - X ) \ds ds 
n 1 .X 


”* 1*2 


Here the unit vector r in direction of line of 
sight has the components (see Fig.l3) r - (cosG, sinG) 


and s(— (s^,s^)) is the vector displacement from the 

• . X ‘ 

instrument sett.ing x. 


/ 


->■ 

At a given value o.E! r (i.e. a given view) S(x,r) (=s (x,r).+ 
noise]! is the recorded image.— A normal two-view stereo image 
would have tv;o images taken at tv70 different values of r* 

/V 

and ^ 2 ) . In the "photographic" mode, the depth of 
field, A^, is typically very large compared to the 
lateral resolution 

F For our example of an asserubly of diffuse targets 


(Eq. 8) , the signal can be evaluated explicitly: 


-1 

n 


f 2 2 \ / 2 2 \ 

1/2 ’ 

)( v -^2 ) 



JXp 






2. A 2 

i "l" A 1 

n 1 


^ 2^. 2 


( 11 ) 

We allov; for measurement error v/xth a noxse term'N(x,r) 
which is to be obtained by constructing independent 
realizations of the viewed field for each viev/ing direction 
r. We assume the noise to be homogeneous (i.e. <^N{x,r)> 

the'. expeatation\= = constant). 

■ Thus the net measured signal is 


,S(x,r) = S^(x,r) + N(x,r) 


(12) 


The Fourier transform of the signal alone can be explicitly . 
evaluated in the form 




a a exp 
n n 


l,:r> , 2,. 2, 

j(k*r) (a^ ) 


n 


X exp (-ik*x ) 
*■ n 


(a^2^A2^) 


( 13 ) 
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The noise field itself v;e take as v/hite, that is. 


= No ® 

k 




U4) 


where is a random quantity uniformly distributed in 
k 

(0,2ir). 

The deduced Fourier transform of the object field 
is obtained by dividing by the transform of the instru- 
mental response function. 


s^^(r) n (r) 
"'k . h * 

*r 


7 ^ = s /g 

k t g^(r) g^(r) 


(15) 


it 


The instrument response function transform is 


=•• exp 
k 


’t '"a 
k'rA^ 


(kxr) A 



(16) 


To avoid singularities at high k we introduce a weighting 
function v/(k) and get our final expression for the processed 
image: 


with 


= 

k 


S-> = 
k 


s_,w_yg 

k k k 


S‘^, (r) + n. (r) 


(17) 


( 18 ) 
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There are two limiting cases 
i) Deconvolution Solution 


This arises if v;(k) 
solution 


1 and noise ->■ 0, giving the 
k 


.# _ 




■ ->■ 


= S Ck,r)/g(k,r) 
k 


ii) Matched Filter Solution 

This arises if w{k) = where w(K)^p^ is the 

optimal weighting function specified by choosing 


■' 

9 

k 


1. k 

’>v) 

k 

g->l 

h| 

> 



(19) 


2 12 
where is the signal power 




'k 


white noise case: 

'2 


-> 

k 


w. 


94. 

k 


2 2 
1 +. N /S 
-- 0 0 


k 


2,2 
+ N /S 
o ' o 


For the 


( 20 ) 


The matched filter solution in the absence of noise n^. 

k 

is then 


^ •>* ^ 

a ' ( k ) = (k , r ) g* (k , r ) 


^1 + (N/S)^) 


(ig (k,r) 


:} + (N/S) 


( 21 ) 
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For the case of viewing directions (21) may be genera- 
lized to a weighted average over the various viev7s: 


a'ih 


- E 


+ (W/S)^) ) 


s^(k,r. ) { ^ 2 

° ^ [ |g(k,r.) 1“ + (N/S)^ 

vxews I ^ ' 3 ^ I 





^ w(k,r . )/ > V7 (k 


( 22 ) 


A fundamental problem in reconstruction of the 
multi-dimensional image arises when there are only a • 

y\ 

discrete set of values of r at which measurements have 
been made. Although it would appear from Eg. (21) that 
a large number of separate viev;s are required to deter- 
mine the spatial function cr(x) in general, experience ■ 
tells us that this : is ^commonly not the case. For^ example, 
the principle that only two two-dimensional projections 
are often all that is required to reconstruct a complete 
three-dimensional image in universally used in conventional 
stereo optical systems. V7e expect, hov/ever, that such ' 
simple reconstruction is possible only under certain 
conditions; that is, when the object is made up of a finite 
number of discrete and identifiable points (or lines in 
3-D) . In order to test this hypothesis vre simulate the 
essential features of such scenes by taking much smaller 
than both A and A^/ and choosing large enough that 

noise can be ignored (n^ -> 0) , 

k ■ 

In this case the measured signal takes the form 


including homo sapiens 



and the Fourier transform of the signal is 
N' 


= Z- % 


1 


a ’ exp - ^ 


n 


x(kx^)2 A,^-ik*J 


n=l 


(24) 


When a sufficiently large number of views have been 
measured it is possible to select r at each k to be per- 
pendicular to ic (r = n(k) where = 0). In this case 

the signal has the form 

n(k))= £ a a ^ exp‘(-k^A ^/4 - iS*x ) (25) 

n=l ^ 


(r = 
k 


This signal in x-space has the form 


s„(5) 


cr(x4-s) G„ (s) ds ds 
2 5C y 


where 


G2(s) = expf-s^/A^^) y/" (■n’A2^) 


(26) 


Thus, when enough views are available, it is possible to re- 
constitute the signal as it v^ould have been seen by -a sensor 
with the orignal lateral resolution A 2 in both range and 
azimuth. This is true regardless of the value of the original 
range resolution A^. This is the essence of the multiview 


process , 



IV. 2, 2 Vector Mode 

• The basic concepts of this section closely follov;s 
that of Section IV. 2.1 v;ith extensions to account for the 
vector field. 


In this mode a velocity vector field is specified 
over the spaticl region: u(x,y) = (u(x,y), v(x,y)). 

For' example calculations^ we have chosen a velocity 
field that simulates a cross-section of the horizontal wind 
field of a storm advancing through the surrounding 
atmosphere (Fig. ,12). This field can be represented for 
simulation purposes by 


u(x,y) 


vortex centers 



i=l 



V (27) 

•> 

where = position of vortex center i 

a^ = 1/e radius of vortex i 
= constant vector velocxty, 

= (0/0,r) and r = vortex strength 
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Fig. 12. Typical Velocity Field Distribution 

For these calculations we take the scattering 
coefficient to be independent of position. 

Eg. 1 again gives the zeroth Doppler moment sensed 
by the instrument at the point x i.e. 


S^(x,r) = G(x^x\r)dx' 


This would be evaluated by the . instrument from the measured 
Doppler spectrum 1(f) : 

~ JlCf)df 

In additon we must sense the first Doppler moment: 


S|| (x,r) = (x' ) u„ (x,r) G(x-x' ,r)dx' 


■; 


fl(f)df 
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( 34 ) 


-y 


S (x,r) -> s (k,r) = (au) ^ g(k,r) 

k 

-y '' -y '' -y ^ 

S^(x,r) ^ s^(k,r) = g(k,r) 


(35) 


The simulated measured data is now tv?o quantities in k space 

-y ^ -y 

at a given viewing direction r^ S^(x,r) and Sjj (x,r) . 


Combining Eg. 32, Eq. 34 and 35 x-je have 

s (k,r) = [(cru)_^ cos9 + (<?v)_^ sinejgCk^r) (36) 

I ' k k 


For this vector field tx^o or more views are required 
to reconstruct the scene. Taking Eq. 36 at two different 

A 

viewing directions r^, and r^ (or at viev/ing angles 6^ and 

0^) \ 7 e can solve for (au) ^ and (av) ^ to give the matched 

k k 

filter solutions for the velocity components as' 


(au) ' 
k 


Sjj(rj^). sine^gg^. 

|gil^ + (N/s)^ 


A ^ 

s^"(.r'2) sine^gg2 

|g,|^ + (N/s)^ 


'sin (0 2-0 


(37) 


iav)\ 

k 


-<s y /N ... * 

-S|j(r^) cos 82 ggj^ cose^gg2 

+ — 


|gj_l^ + (N/S)^ 


g.l^ + (N/S)^ 


js±n{Q^-Q^) 


(38) 
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where g = 1 + (N/S) . These components cen be averaged 

over a number of sets of pairs of give 

CEq, 37) (39) 


(Eq. 38) (40) 

ft 

Note that because of possibility of small denominators in 
Eq. 37 and 38, these sums should be restricted to only the 
more widely separated pairs of views. As before logarithmic 
interpolation can be used to interpolate in regions of 

small (au) ' and small (av) ' . Performing the inverse transform 

. k 

to X space gives • 

u(x) = (au)'/o'' (41) 

and 

v(x) = (0v)//a' (42) 


k 


I 


... 1 

N 

pairs . - 

^ pairs of 

views' 


Cav>4 

k 


N 


pairs 


E 


pairs or 
views 
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IV. 3 Simulation of a Multiview System 
IV. 3.1 Introduction 

The procedures described in the preceding section 
have been formulated on a computer and a number of simula- 
tions carried out for both scalar and vector fields. The 
basic features of the output and display are shown in an 
example in Fig. 14, 15. Here the scene is taken to be 
the v7ord NASA represented by a sequence of simulated 
Gaussian shaped scattering centers. Tv;o types of output 
display are shown: a contour plot of the processed i ^age 

field and a raster plot where the intensity is proportional 
to the computed local image density. In Figures 14a and 15a 

the original scene is shown and in Figures 14b and 15b the 

\ 

scene as seen by a single ^^ensor v/ith a range resolution 

i 

equal to 1.^4 times the image letter heights. Matched filter- 

1 

ing has been used to construct the images v;ith an assumed 

I 

SNR for the matched filter of 1.0. Hov;ever, no noise vias 
actually present in the simulated data. For a single viev7 

i 

(Fig, 14b, 15b) the sensor resolution is inadequate to 

i 

1 

resolve the pattern. Figures 14c, d, e and 15c, d, e sho;>7 

1' . ■ ■ ■ 

how the processed spatial resolution improves and the 

. i ' - 

pattern emerges as the numberl of views is increased. The 

i 

pattern is recognizable with jfour viev?s and very clear with 

I 

24 views. 1 
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Example of a Multiviewed Scene - Contour Plot 


a) 


b) 


c) 


d) 



1 View 




1 \ '^^\\ \'\ \ 

.xM|^ 

> ^Wc 7 

^ \ ' 0^5us/>a > 

^ w . 

fedcfi^Sr 

37 pC^ 

ja ^ 'V^- 

r^w'A r \ 

^ X/. ^ . 


Figure 15 







The b^sic problein nov7 is to quantify the processing 
procedures so that the relationship between single sensor 
resolution, number of views, processing filters and final 
image resolution can be established. A quantitative mea- 
sure of the quality of the image reconstruction is required. 
The Rayleigh criteria for defining the resolving power of 
an optical system is a standard method. It gives a good 
measure of the smallest scale structure that the instrument 
v/ill sense. However, it is not very useful in this applica- 
tion since it gives no indication of distorted or false 
images or-'ofthe presence of other processing induced irregu- 
larities. V7e have found the cross correlation coefficient 
between the reconstructed image and the original scene to be 
a useful measure of image quality. This quantity is defined 
as 


C = 
o 


Jo'ix) a^Cx) dxdy 


^c^(x) dxdy ya^(x)dxdy 


% 


V7here a^(x) is the original scene and (x) the reconstruction 
For simple localized images, this correlation coeffi- 
cient can be interpreted directly in terms of a spatial 
resolution. For exainple, for a source distribution of the 

^^2/ 2 -v^// 2 

•form e '^^s and an image distribution of the form e 
the cross correlation is 


2 


(44) 



For an image ^^7hich is much broader than the original 

2 

source the quantity 2/C^ is effectively the ratio of the 
image dimension to the original source size (A^/A^) . 

Because of this simple relationship between resolution and 
correlation coefficient we will use the quantity as 
the appropriate measure for comparing different processing 
techniques. 

IV. 3, 2 Results 

To test the multiview process a number of simulations 

ware run using test patterns for both scalar and vector 

fields for the original scene (see Sec. IV. 2). The data 

is presented in two forms; contour plots and raster plots, 

the former being more useful for quantitative analysis 

and the latter giving a more pictorial sense of the data. 

In Figures 16 and 17 v/e shov; the effect of varying the 

number of viev7s for tv/o different sensor resolutions for 

a scalar field. The original scene of three Gaussian 

shaped scalar scattering sources is also shown. The upper 

four figures correspond to a very poor initial spatial 

resolution (2A^ = AX = 24 units) and show a substantial 
1 o 

improvement of the scene resolution as the number of views 
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1 view 


4 View 


2 View 


8 View 


Original 

Scene 


2 View 

2A, = 8 


1 View 
2A, = 8 


Fig. 16. Reconstruction of a Scalar Field (Contour Plots) 


















is increased. The lower three figures have a higher spatial 

resolution to start with and achieve the same resolution 

as the first set at a lower number of views. This effect 

is shown more quantitatively in Figures 18 and 24. In 

2 

Figure IS we show how the correlation quantity (for a 
two-view model) decreases with deteriorating sensor reso- 
lution, and shows roughly the inverse relation expected 
from Eq. 44. In Figure 24, we show in a quantitative way 
how increasing the number of viev/s can offset this degrad- 
ation. 

Two signal-to-noise values enter into these simu- 
lations: a) the actual instrument noise that is to be 

added to the simulated measured quantity ‘ before processing 
and b) the assumed or expected value of the signal-to-noise 
value that is to be included in the matched filter expression 
In most of the calculations shown a rather low {but finite) 
value of instrumental noise v^as selected in order to demon- 
strate the expected profile of the signal. Various values 
for S/N of the matched filter were selected for demonstra- 
tion purposes, ^o^vever, it is to be noted that in practice 
the best results will be obtained by matching the filter 
S/N to that of the actual signal- In Figures 3.9 and 20 
we s-hov7 the effect of varying the actua.l instrumental or 
background noise in a typical example. Here it is seen 
that too lov; a value of S/N in the filter leads to a 




Ai 


2 

Fig. 18. Image quality parameter Cq for a tv/o-vieW’ 

■ reconstruction of the scalar field shown 

in Figs. 16 and 17- = 1/2 range resolution 
betv/een 3 db points) . Values for the (voltage) 
S/N level assumed in the matched filter are 
indicated. 
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ASSUMED SNR OF MATCHED FILTER (DB) 


Figure 20. Variation of the image quality with 
assumed value of the signal-to-noise 
ratio in the matched filter. (Pictures 
in Figure 19 refer to. points at a, b, 
c , d and e) 
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broadened pattern v/hereas too high a value leads very 

noisy looking data. The 'best' image, as defined by the 

2 

correlation coefficient being maximum, occurs when 
the S/N of the filter is matched to that of the signal 

(Figure 20) . 

In Figures 21, 22 we shov; a similar set of simula- 
tions to those of Figs. 16, and 17 except nov/ for a vector 

field. The scalar quantity presented in these figures is 

2 2 1/2 

now the speed (u 4-v ) ' rather than the cross-section. 

The vector field chosen for this example is that given in 
Section IV. 2. 2, Eq. 27 and v;as chosen as a rough simulation 
of the horizontal velocity field of a moving storm. 

In addition to the contour and raster plots of the 
speed, we also show vector plots of the velocity (Figure 23) 
The velocity vector at any point is given by a line whose 
length is proportional to the speed and whose direction is 
indicated by 'the point at the front of the line. 

From the tv/o-diraensional reconstructions just des- 
cribed we have constructed a set of curves relating the 
number of viev/ing directions to the basic sensor resolution 
length (2A^) for fixed values of the correlation coefficient 
C^. The plots are presented as a function of the inverse 
of since this is dire'r’tly proportional to a quantity of 
direct interest to the electronics design - the peak data 
handling rate. 




24 View 

A — 0 A 


8 View 


4 View 


2 View 


8 View 


2 Viev; 


4 View 


Original 

Scene 


Figure 21. Contour plots of image reconstructions for a vector field. Quantity 
plotted is the local magnitude of the velocity. (All Graphs S/N = 1, 
HI-PASS FILTER) 














Figure 22. Raster plots of image reconstructions for a vector field. Quantity plotted is 
■ local ■ gnitude of th._ velocity. (All Graphs S/N = 1, HI-PASS FILTER) 




























The qualitative features of the images shov/n in 
Figures 16. to 23 appear generally consistent with the 
quantitative correlations shown in Figures 24 to 26, 

For the scalar mode the number of views must be increased 
roughly in proportion to an increase in resolution length 
to maintain the same image resolution. Of course, a sat- 
uration in this relation is reached v/hen the reconstructed 
image is close to the original scene. 

The vector field reconstructions (Figure 21, 22 and 23) 
are not as satisfactory as the scalar field reconstruc- 
tions. The same general trends are evident, but substan- 
- tially larger numbers of views appear to be required to 
fully resolve the vector image. It is felt that a more 
detailed analysis of the filtering and reconstruction pro- 
cess for the vector field v/ill be required to develop an 
improved algorithm.. Preliminary attempts to evaluate the 
efficacy of various filtering techniques to improve the 
image quality are discussed in the next section. 
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NUMBER OF VIEWS 





Figure 2 4 . Number of viev7S' required to maintain a given 
processed image resolution as a function of 
the single sensor range resolution (3 db 
resolution = 2A) . The parameter is the 
■ cross correlcition between the processed scene 
and the actual scene. (Scalar Mode) 


NUMBER OF VIEWS 



Figure 25. Number of viev;s required to maintain a given 
processed image resolution as a function of 
the siTig]e sensor r£inge resolution (3 db 
resolution -- 2A) . The parameter Cq is the 
cross correlation between the processed scene 
and the actual scene. (Scalar Mode) 
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NUMBER OF VIEWS 


(S/N= 1 . 0 ) 



Number of visv/s required to maintain a given 
Figure 26 . processed image resolution as a function of 
the single sensor range resolution (3 db 
resolution - 2A) , The parameter is the 
cross correlation between the processed scene 
and the actual scene. (Vt^ctor Mode) 



IV. 3, 3 Image Enhancement Procedures 

The essential feature of the multiview process v?hen 
the number of viev/s is finite but at the same time the 
sensor resolution is poor is that the measurement incom- 
pletely samples the scene. In some directions detailed 
information is available but in others very little. In 
the Fourier (v^avenumber ; (k ,k )) domain this means that 

^ y 

data is given only along certain rays emanating from the 
origin (k^ = k^ = 0) . A variety of methods for filling 
in or interpolating this data can be discussed and v;e 
v;ill present one preliminary approach to this problem 
later in this section. First, however, v;e consider the 
proper choice of filters v/hen the sampling is adequate. 
Simply adding together the independent viev?s vzill not 
yield a perfect scene reconstruction. A point function 
(a (x) = 6 (x) ) v;ill yield an image o' (x) in this process 
which decays inversely as the radius {o' ~ 


2 , 

X -1-y ) . 


A very simple filtering procedure has been used in 
Figs 27/29 to improve the image resolution of such data. 
This is basically a high pass filter. The final pictures 
are less smooth but the spurious x-/ide wings associated 
with the multiviev; superposition are eliminated. 

To investigate how v;ell filtering procedures could 
improve the- I’ccons truction of an incompletely sampled 
scene a single Gaussian source was examined (see Figure 30) 
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Ficure 27. Effect of high pass filtering (scalar mode, all graphs N/S =1). Note 3rd and 
<lth column graphs are Raster Plot representations or 1st and 2nd columns. 



















Figure 28, Effect of high pass filter (vector mode, speed distribution; all graphs N/S = 1) 
Note: 3rd and 4th column graphs are Raster plot representations of 1st and 2nd 
columns. 























An extreme case v;as taken with a full v;idth instru- 
ment longitudinal resolution of = 24 units and two views 
The instrument noise level was initially taken as zero 
(Fig 30a and 30b) and then taken as 14 db (Fig- 30c and 
30d) . The 2nd row of pictures shows the reconstructed 
scene after applying the matched filter. This would be the 
same situation as for all previous pictures (except with 
the addition of the application of the hi-pass filter) . 

The 3rd row of pictures shows the scene after a 
post processing (described below) and finally after a 
hi-pass filter was applied (4th row of pictures) . 

•Essentially the post-processing procedure attempts 
to fill in regions in k-space away from- viewing directions 
where the da-ta is poorly represented or non-existent. 

First the data in k-space (i.e. the actual data as 
measured from the instrument) is homomorphically filtered, 
(i.e. the log of the data is taken and the imaginary part 
(the phase of the. original data) is smoothed. Each point 


in k-space is then weighted by the quantity 


where. W = / ’ g? (1 + N^/S^) / (g^ t ) (see Eq. 20) . This 

O X 

vxsws * 

factor is a function of viewing angle, and correctly retains the 

integrity of the data along viewing directions. In 

addition to giving a unit \- 7 eighting to directions near 

the viewing directions it assigns low weighting v/here the 

da-ta strength is very low* : The data is then averaged in the 


y ^ o 

^ views 


angular direction according to these v/eighting factors, 
i.e. unit vreighted ('good*) data are not changed, v;hereas 
data with low v/eight is interpolated from more highly 
weighted data. This has the effect of filling in data 
in poorly sampled regions from data in well sampled regions 
This procedure works reasonably well for the single 
Gaussian case (Figure 30). The same procedure was applied 
to the standard three Gaussian source case. For the case 
considered of tv;o viewing directions and 2A^ = 24 it was 
fo^Ind that the internal structure of the scene xvas not 
very vrell represented although the procedure did strongly 
localize the image structure and suppressed the spurious 
far distant contributions. (Figure 31) 

Based on these preliminary results’ it appears tha-t 
substantial image improvement can be achieved by the post 
processing procedures. Hox^ever, the treatments reported 
here are not universally applicable to different scene 
types and more experimentation is needed. 



Reconstructions after applying matched filter 




Reconstructions after post-processing 
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V. Conclusions and Recommendations 

The essential conclusion of this study is that 
aircraft flight velocities are high enough and severe 
storm lifetimes are long enough that it is possible 
for a single airborne Doppler system, operating at a 
range of less than about 20 km, to view the storm area 
from two or more substantially different aspects before 
the storm characteristics change appreciably. This per- 
mits resolution of the spatial distribution of the velo- 
city vector by a single coaxial system operating in a 
fly-by mode. In addition, it is suggested that, when 
the instrumental range resolution is poor, post- flight 
processing of such fly-by multi-view data can be carried 
out to materially improve the spatial resolution of. the 
data. 

When peak laser power levels and/or peak data 
handling rates are limiting problems (as is expected to 
be the case), substantial gains in signal-to-noise as 
well as velocity resolution are possible by implementing 
the •multiview mode V7ith strongly degraded single sensor 
resolution. Here the . long pulse times permit a higher 
ratio of coherent to incoherent integration and a more 
accurate determination of velocity. The penalty required 
is that of a more complex optical scanning system as v7ell 
as additional real-time or post-flight time data storage 
and processing. 
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Tv 70 essentially different scanning configurations 
have been identified, one a long-range (out to 10-20 km) 
roughly horizontal plane mode intended to allow probing 
of the velocity field around the storm at the higher 
altitudes (4-10 km) . The other is a shorter range (out 
to 1-3 km) mode in v;hich a vertical or horizontal plane 
is scanned for velocity (and possibly turbulence) , and 
is intended for diagnostics of the lovzer altitude region’ 
belov; the storm and in the out-flow region. The first 

mode has a prototype form and an advanced form. In the 
prototype, only a single horizontal plane v/ould be scanned 
for the u and v components of the horizontal velocity as 
well as certain combinations of their gradients (vorticity, 
convergence) . In the advanced mode, multiple plane (simul- 
taneous) capability is desired (10 or more stacked hori- 
zontal planes) and v/ould be used to extract u, v, and w 
velocity information (and gradient) by processing the 
measured u and, v data. 

Application of the advanced mode may be primarily 
to complement the dual-Doppler radar measurements because 
of the limitations imposed by cloud attenuation as well 
as low-altitude attenuation. These may prevent the laser 
system by itself from acquiring all the data required in 
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the full storm-related volume to permit reconstruction 
of the full velocity field. 

It is recommended that detailed design considera- 
tion be given only to the short-range mode and the proto- 
type long-range mode. Implementation of the advanced 
mode should be delayed until successful testing of the 
prototype mode is accomplished and/or additional analysis 
of the full 3D mode has been carried out. 

Both the short-range and the long-range modes 
will require multiple-viev; capability.; that is, scanning 

of the viewed area from two or more different aspect 
angles (within a possible range from 10® to 90°). These 
multiple views are required for two purposes: for ex- 

tracting velocity vector data and for enhancing range 
resolution. A system that viev/s each point in space from 
one direction only can only be used in conjunction with 
some assumptions such' as horizontal uniformity of the 
wind field (as in VAD systems) and is not one of the 
permissible modes for this application. 

Tradeoffs between number of views and range reso- 
lution exist. A degraded range resolution (due to a 

pulse time) in the real-time data may be ' acoeptable , 
if a sufficient number of different views exist to allow 
post-flight processing for image reconstruction. Detailed 
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plots of this tradeoff are given in Figure 24, 2S, and 26. 
These data are derived from computer simulations and give 
the interdependence of number of views and sensor resolu- 
tion that must obeyed in order to maintain a given processed 
image quality. • 

Under circumstances where pulse lengths are short 
enough to produce adequate range and velocity resolution 
without processing, only tv?o views are required for the 
vector velocity field. However, it must be remembered 
that the vector fields are basically being differenced 
in the processing at least once and perhaps twice (once 

to get the components and again to get the vorticity 
and/or • convergence) . Thus high signal-to-noise levels. . 
are required as well as a finer grid. 

Even for a many-view system the averaige output 
data rate (the data required to define the flow field 
that must be recorded per second) is relatively low 
(less than 100 data points per second) . Without some 
real time processing, a 10 to 20 times higher recording 
rate may be required for the optimal multi view, continuous 
scan system. The advanced system V7ill also require a 
lO-to-40-times higher rate. Hovrever, even here the total 
amount of data to be stored during a single measurement is 
expected to be less than 3 x 10^ data. Thus the data 
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storage is probably not severe. However, the peak dat4 

f 

handling rate at the front end may for some designs be 
very high with a very low duty cycle. Thus it is recommended 
that serious consideration be given to modes in which the 
peak data handling rates are traded off against fairly 
extensive post-flight processing procedures. 

Design Recommendations 

It is suggested that aspects of the following 
problem areas be addressed in the design analysis. 

(1) Registration and storm region transport . Combining 
multiple views requires that the data be recorded in a 
fixed spatial reference frame rather than just the a/c 
frame. Since different views are displaced in time, the 
relative displacement of the a/c to the storm needs to 

be accounted for in the recording of the data. 

(2) Techniques for obtaining the zeroth, first, and 
second Doppler moments (backscatter level, mean velocity , 
and velocity spread) in each space-time resolution element 
from the Doppler spectra.. 

(3) Scanning t echniques for obtaining 2 , 3, 4 and/or 
many (10-3-0) •aspect looks at the target region in a 

fly-by, 

100 



(4) Relation betv?een peak front-end data handling 
rates and pulse duration and shape. 

(5) Tradeoff between pulse duration and number of views 
(i.e., betvreen peak front-end data handling rates »and 
multi-view scanning techniques) . 

(6) Real-time partial processing and display for purpose 
of data quality monitoring. .. and reduction of data storage. 
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MATCHED FILTER PROCESSING - Ali EXAMPLE 

In order to demonstrate the nature of the matched 
filter process, v;e consider 'the detection and processing 
of the signal resulting from the illumination of a point 
scatterer. For convenience we- treat only the one-dimensional 
problem. 

The scatterer is located at the point x = x^ and is 
described by the distribution function 


h (x) = aQ6(>:-XQ) 


( 1 ) 


For this example we choose a Gaussian shape for the instru- 
ment response function: 


G(u) = exp(- u^/A^)/f /rll) 


( 2 ) 


The detected signal is 

S (x) = ^CT (x' j G ix-x ' ) dx^ -f noise 
and its Fourier transform is 


( 3 ) 


S (k) = a (k) g (k) + n (k) 


( 4 ) 


For reference v;e note that the transform of the original 
scattering function is 


a { k ) = a ^ exp ( i kx^ ) 


and that of the instrument response function is 


( 5 ) 


In the total absence of noises the signal may be 
prot.essed by a deconvolution filter to yield the original 
scatterer distribution: 

ilex 

cr'(k) = [sCk)J = >(k) = o^e ° ( 7 ) 

In the presence of finite noise however, this 
deconvolved signal becomes 

0' (k) = 0 (k) + n (k)/g(k) ( 8 ) 

In regions where g(k) is small (k >> 1/A) the noise is 
greatly amplified and can effectively destroy the imaging 
process. Here an appropriate matched filter should be 
used. For optimal processing in the presence of white 
noise, the deconvolution filter (g(k)) ^should be multi- 
plied -by a weighting factor V'?(k) given by . ' 


w(k) 


|gl^(l+(N/S)^) 
(N/S)^ + !g|^ 


(9). 


where S/N is the expected signal/noise ratio 
Thus the processed signal is 


w _ s (k) q* (k) (1-KN/S) ^) 

|gj'^ -P CN/S)/^ 


(10) 


and the noise is 


n' (k) = 


: n Cklg^ Ck) (1 4- (N/S) 

\g\^ + (N/s)^ 


( 11 ) 


For our point scatterer example this signal' has the explicit 
form 


a' (k) = c^e 


' ' (1 +' '(N/S) • • 

1 + (tys)^ exp(k^A^/2)) 


(12) 


In x-space this processed signal has the form 


o (x) 


a 


o 


ik(x-x ) ' 

e - • (1 + (N/S> ) dk 

^ + (N/S) ^ exp(k^A ^/2)j 


(13) 


For high signal-noise ratios (N/S << 1) an approximate 
evaluation of this integral can be obtained by setting the 
denominator equal to 1 for 1 kA j < /^n (S/N) and to infinity 
otherwise. In this approximation 


4or 


sxn 


.r X-X 

r o 


2 A*. 


(x) = ^ 


TT j 


2 j 



(14) 


wnere 


A* = 


TT, 


A 


^ /^n~TWr 

Thus the width of the processed image is reduced from the 
instrumental width by th e factor 


4 /in 


The extent of the image improvement that can be 
achieved at high signal/noise ratios depends explicity 
on the shape of the instrument function. For example, 
a box function of the form: 


G (u) - 1 for [ u| < A 


= 0 otherwise 


has a reconstructed image of the form 


G '(x) = a: 


ik(x--x ) „ 

f • (1 t (N/S) dk 


1 + (N/S) 


2/ kA 
i sin kA 


An approximate estimate of this integral (as before) leads 
to a simple but different dependence of the processed image 
width on the (voltage) signal to noise ratio: 


n/i/2 (S/N) ) = tt/^T(s7nJ~, 

n ■ yOvv 


pov/er 


This box shaped resolution function exhibits a much greater 
reduction in image v;idth at high S/N values than does the 
Gaussian form. This results from the sharp cutoff at the 
edge of the resolution function and the resulting relatively 
rich wavenutnber spectrum. 
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In Figure 32 v;e show numerical calculations of 

I ■ ■ 

Eq. 13 for the expected profile of the matched filter pro- 

I : ■ ■ ■ 

cessed signal for various values of signal to noise using 
the Gaussian response function (Eq. 6 ) • The figure shows 
the rapid narrowing of the resolution function as S/N . is 
increased at low S/N values changing to a much slo;^er 
dependence at high S/N as would be expected from the 
logarithmic dependence given in Eq. 18. These profiles- 
are the 'expected values' for the signature. Any given 

measurement with finite noise will appear more irregular. 

% 

In general# use of matched filters %-7hose S/N values are 
less than those of the actual signal leads to smooth 
looking output with a clearly defined but relatively broad 
profile. However# the optimum matched filter (for which 
the value of S/N in the filter is ma-tched to that of the 
data) always provides the best localiza-tion of the source. 



SIsm =. 20 db 
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r^atched filter reconstructions of, the profile of a point target at 
different signal to noise ratios. The sensor is assumed to have a 
Gaussian response function in range. The 3 db resolution widths 
(relative, to the width, expected in the small signal limit) are 
indicated . ‘ . 




APPENDIX II 

CLASSIFICATION AND DESCRIPTION OF WIND FIELD 

CHARACTERISTICS OF SEVERE LOCAL STORMS 

I. INTRODUCTION AND SUI#IARY 

One could classify severe storms according to the 
overall size of the phenomena v/hich produce severe effects 
into severe local storms (thunderstorms or tornados) and 
into severe extended storms (hurricanes or high pressure- 
gradient storms). However, in the following the notion 
of severe local storm , is applied to dangerous and 
damaging phenomena occurring in all of these storms because 
the severe features of such storms are usually v;ell localized 
in space and time. Actually, the notion of a .severe local 
storm has a restrictive meaning for the meteorologist and 
weather I forecaster . .when lightning strikes practically 
continuously within and between the clouds of a hugh 
thunderstorm..." (ref.l) or, more definitively, a so-called 
"supercell thunderstorm" v/hich, by definition, produces hail- 
stones of the size of tennis balls in temperate latitudes 
(p. 159, ref. 2 ) . Dangerous situations connected v/ith the 
motion of air around and in storms are caused primarily by 
relatively sudden and large changes in the velocity and in 
the direction of local wind features . Such changes may be 
related to: (a) strong down-drafts and up-drafts around and 

in thunderstorms or in large gusts / (b) excessive wind shears 
from strong gusts in pressure-gradient storms, from tornados. 


and from cold outflows of thunderstorms, and (c) obstacles 
in the path of strong winds. In this Appendix, synopses are 
for the categories of storms and of features around and in 
storms v/hich produce such phenomena. A measurement system : 
must be designed to conform to the space-time velocity 
structure of these storms. An overall summary of these 
dimensions is given in the following section. This summary 
is based on a more in depth discussion that follows. 


I.l Overview of the Space-Tiire Velocity structure of 
severe local storms 

210 Thunderstorms (air-mass, single-cell) 

211 Depression / trough and orographic thunderstorms 

life cycle; 1/2 - 3 hrs. 
cell“size; length: i5 - 50 km 

width: 5 - 15 km > . 

height: up to 12,000 m from 400-500m 

above ground _ • 

approach speed: 15-60 km./hr (stationary at mountain 

crests) 

severe gusts: 15-25 ra/sec : cold-outf lov7, roll-cloud- 

gusts: 5-10 m/sec increasing up to 8-10,000 m 

severe updrafts; 5-20 m/sec forward cloud portion 
severe downdrafts:' 10-40 m/sec backv/ard lo\'7er cloud 

position . • 

O 

strongest shifts: near ground at front (up to 180 ); 

. throughout storm, might veer ; 

up to 90 immediately after storm 
strongest winds: 20-40 m/sec at 300 to 500 m near ' . 

hail area; can double at mountain crests 
212. Multi-cell and area thunderstorms 

life cycle: 1-5 hrs. , 

cell-size: height: 50-rl50 km 

v/idth: 100-30 0 kra 
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approach speed: is-60 km/hr (stationary at ' 

mountain crests) 

severe gusts: 30-60 m/sec and largest in major cold. 

outflov/s, particularly v;hen supported 
by down slopes; at front-base and • ' 

l0v?er front of roll clouds. 

severe updrafts: 10-25 m/sec in for^vard cloud portions 

severe downdrafts: 15-50 m/sec in backv/ard lov/er 

cloud portions, enhanced by down 
slopes. 

strongest shifts: same as under 211. . 

strongest winds: 25-50 m/sec in cold outflovzs 

213. Supercell thunderstorms • ' 

life cycle: 1-10 hrs. ' . 

cell size: length; 50-150 km . 

vjidth: 20-50 hrs * 

approach speed: same as under 211 

other features: same as under 211. 

214, Thunderstorm squall lines (prefrontal) 

life cycle: 1 -10 hrs 

cell size: length: 20--60 km (in direction of motion) 

v/idth: 250-500 kra (length of band) 
approach' speed: 15-60 km/hr • 

severe gusts : 20-30 m/sec at front of cold-air 

outflov7S under leading edge (roll cloud) 
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severe updrafts and downdrafts; simiXar a‘s or more 
intensive than 211. 

shifts and strongest v/inds: similar as or more iritensive 

than 211. • 

22. Tornadoes . 

typical structure;, interconnect clouds and ground; .• 

electrical 'discharge. • 

path-length; 0.5-15 km (150 to 500 km 

on rare occasions) ' . " 

path-v/idth; 10-300 m (0.5 to 10 km on ' 
rare occasions) 

approach speed: 15-30 km/hr (50 to 80 km/hr on rare 

• ' ' . . < . .... 

* ‘ occasions) . * 

v/ind spaed: 20-70 m/sec (100 to 200 m/sec on rare 

occasions) 

updrafts; suction vortices 30 to 100 m/sec 
persistance: 20 to 50 min in families* • 

recurrence: 30 to 60 -.min intervals . , 

23 . Funnel Clouds (cold air) 

• no tornado potential _ ^ • 

typical structure (connected to clouds^ may not 
touch ground) : V 

path-length: 1-4 km (up to 10 km on 

• rare occasions) 
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path-width; 1-5 m (up to 20 m on rare 
occasions) 

approach speed; 15-50 km/hr 

v;ind speeds: not really knov/n; :: 2-10 m/sec 

persistence: 3-10 min (15-20 mins on rare occasions) 

Dust devils . 

typical structure (rising from the ground; no dis- 
charge) 

path-length: 0-5 km (up to 10 km on 

rare, occasions) 

path-v/idth: 10-100 m (up to 500 m on 

rare occasions) 

heights; 30-200 m (up to 1000 m on rare 
occasions' V7ith disturbances to 
2000 m) 

approach speed: 0-20 km/hr (raove with average wind 

if > 5 km/hr) 

v;ind speeds: not really knov/n; n 5-30 m/sec (up to 

50-100 m/sec on very rare occasibns) 
updrafts: not really knovm; 10-50. m/sec (up to 

100 m/sec on rare occasions) 
persistence: 2 to 5 min (30 mins to 4 hrs on rare 

occasions) 

Density c urre n t vortices 

Powder snow avalanches 

bulk-flow velocity: 30-100 m/sec 


vortex-f lov; velocity: ~ 2 x bulk velocity (up to 

200 in/sec) 

vortex duration: • 2-3 sec for each cell 

vortex cell size: ~ 0.5-30 m (estimate) 

252 . Nuee ardente 

bulk-flov7 velocity: 3-10 m/sec 

vortex-flow: 70-100 m/sec 

vortex-flow duration: not really known (~l-20 sec 

, estimate) 

size: not really known (0.5—30 m estxmaL.e) 
structure: many short (O. 5-30 metres minimum) hori- 

zontal vortex cylinders. 

253. Sand and dust storms (haboob) 

bulk-flow velocity: 5-20 m/sec (and larger on 

occasion) 

■ vortex-flow velocity: 20-100 m/sec (and larger on 

occasion) 

vortex— flow duration: not really known (~5 to 60 

• min estimate for tornado-like 

structures) 

vortex cell-size: 10-300 m (and larger on occasion) 

254. Typical vortex structures in density currents 

2541. Surge heads 

2542. Rising tunnels 
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2543 

2544 

3. 


31. 

32. 

33. 


Vortex pairs 
Vortices from clouds 

Typical strong V7ind shifts and gusts 


V ^ - (1.44 ± 0.11) V . . 

gust ' wind 


{ 


horizontal 

vertical: above ~ 30 m 


persistance: (horizontal) 1/2-30 sec 

(vertical) ~ 1/3 x (horizontal) 
distance for increase to maximum velocity: 15-70 m 

wind shears' 1/10 - 3 sec ^ 


directional drifts: 15-30 ° near ground 

(up to 160“ observed 'under 
roll-clouds) 

Wavestreaming and rotor streaming 

shifts: 45-160° particularly under rotor- clouds 

velocity of gusts : depends on average v;ind ( -V's, 

above under 3). 

severe downdraf ts: belov/. and near front of roll 

clouds and near lee of obstructions 
Turbulence in the lee of mountains 
(qualitative features) 

C'irc^letion gust cells 

wind shifts (northern hemisphere) : veered with 

increase, backed vrith lull/ 15-30° (up 
to 45° on rare occasions) 
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ceXL sizer 100-2000 la for small to medium size 
eddies 

persistence; 5-50 sec 
extension; up to 1000-2000 m 



2‘. TYPICAL SEVERE LOCAL DENSITY STORMS 

Under this heading, phenomena are combined v/hich 
occur primarily in conjunction with instabilities caused 
by the mixing of air-masses with different temperatures, 
densities and moisture contents. Thus, thunderstorms 
and their associated features, including tornados and funnel 
clouds, are considered; but also dust-devils, v/ater spouts, 
airborne avalanches of pov7der snov; and of volcanic m.aterials 
(nuees ardentes) , and sand and dust storm.s (hahoobs) belong 
here. Cold winds down large mountain slopes (katabatic 
v/inds) , like cold foens, cold Santa Anas, or boras, cause 
similar effects as cold air outflovzs from large area thunder 
storms and are included under vortex-generating density 
flov;s. V7 arm katabatic V7inds (vjarm foens or vjarm Santa Anas) 
are not incorporated here since their turbulent features 
(associated primarily with mixing of upper-level, high- 
velocity winds with ground-level air masses) are usually 
caused by topographic obstructions. 

• ( 2 - 7 ) 

2,1, Thunderstorms (A^ir-Mass, Single-Cell) 

Air-mass thunderstorms (single cell variety) will 
be ignored.. They grow among other large cumulus and 
cumulonimbus in very unstable airstreams (thermals) and 
are usually not severe . 
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(2-7) 

2.1.1, Depression^ Trough ^ and Orographic Thunders torins ^ 

Occurrence : Areas of depressions and troughs v/here 

usually strong tendencies for ascent of moist air masses 

I ; I . 

exist. Ascent of v;arm and unstable moist air masses at 

mountains can lead to (orographic) thunderstorms. Outbrea^ts 

, *1' 

are widespread (especially in summer) . One storm cell tends 
to build another on the back of the very cold air v/hich 
cascades out of the parent storm and spreads around it. 

Groups of cell form^ each nev7 cell taking over, v/here a 
previous cell dies out. The chains of cells can coalesce 
into area storms which may cover areas of- 15 to 30 km across 
and can result in thunderstorm highs . - . 

Typical Cell Structure: (See Pigs. 2.1a and 2.1b.) 

length: 15 - 50 km 

width : 5-15 km 

life cycle: 1/2 - 3 hrs 

speed of approach: 15' - 60 km/hr; center often stationary 

near mountain crests. 

severe gusts.: 15 - 25 m/sec at the front of cold-air outflow 

and near the base in front of the cloud: roll 

cloud. 

severe updrafts : severeness of gusts increases to 5-10 m/sec 

v/ith altitudes up to 8 - 10,000 m. , 

5 - 20 m/sec in forv^ard cloud portion. 

severe dov/ndraft: 10 - 40 m/sec in backward lov/er cloud 

portion (hail area) . 

strongest wind shifts : At and near ground at front of cold 

outflov7 veered by almost 180° . Wind 
veers throughout storm; sometimes up 
to 90° to direction of approa-Gh 
immediately after storra passage. 


strongest winds: 20 - 40 in/sec flovjing out from and near 

hail area at altitudes from 300 to 500 m. 
Wind velocities can double above mountain 
, crests. 

Special Features : Sudden outbreaks of large thunderstorm 

areas can result in thunderstorm highs due to replacement 

I ” ' 

of v;arm air Ijy cold air. These highs may set up pressure 
surges causing sharp and large directional v;ind changes 
(often backing up to 180*^) , as v/ell as violent gusts from 

; 'I 

2 to 15 m/sec many miles away. The location of these changes 

I 

depends on the movement of the pressure surges (at and 
near ground level) from the thunderstom high. Cold outflows 
may affect areas up to 100 av;ay from large thunderstorm 
areas (precipitation dovmdrafts) , particularly if mountain 
terrain supports the descent of the cold masses. 

For a model of hail-producing thunderstorm/ sea Figure 2,2. 



h 


10,000 




h/m 


Anvil Dome 


Rear Anvil 


1 * * • •• • • 

. / • • • * . . • ^ . • * • \ iW 



Froixt Anvil 


mall Hail 



Snow 



Roll Cloud 


. •* . . 




'10 km 


20 km 


Rain 


Large Hail 


Squall Front 




Streamlines Relative to Storm 
Hailstone Trajectories . 

Condensation in Updraft 

Precipitation which has Left Updraft 


Figure 2.2. Typical Model of a Hail-Producing Thunderstorm. 
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2.1.2. ■ Multi-Cell arid Area Thunderstorms ^ ^ ■ 

Occ^tre'nce : Tropical cyclones and reinforcement 

of remains of tropical cyclones; touchdov^ns of jet stream 
at troughs. Strong convective updrafts of large and unstable 
moist air masses at large mountain ranges can result in 
(orographic)’ area and multi-cell thunderstorms. 

■ ■' Typical Structure : May depend on terrain, 

length: ‘ .50 .150 kra 
vzidth: 100, - 300 km' 

life cycle: 1-5 hrs ’ 

spaed of approach: 15 - 60 krn/hr; storm centers often 

stationary at mountain crests. 

Severe gusts: 30 - 60 m/sec and larger in major cold ■ 

outflovjs, pEirticularly v?hen supported by ' 
dov7n slopes of terrain. Also at front- 
base and front of clouds (roll clouds). 

severe updrafts: 10 - 25 m/sec in forv/ard cloud portions. 

severe dovmdraf ts : 15 - 50 m/sec in backward lower cloud 

portions and hail area; enhanced by 
down slopes of terrain (katabatic winds) . 

strongest v/ind shifts: At and near ground level at the 

front of individual cold outflov?s 
and at the rear of individual cells. 

strongest v?inds: 25 - 50 m/sec in cold outflows . 

Special Features : Usually many and strong hail areas v;ith 

hail stones from 3/4“ to 1“ size; often tornados during and 
funnel clouds during and after passage of severe multi-cell 
and are a storms - Development of thunder stG.rm-high enhanced 
emd effects mentioned lincler 2.1.1 apply to ci larger degree. 


2 . 1 . 3 , Supercell Thunder storms 

Occurrence : Same as for multi-cell and area thunder- 

storms V7hen all circumstances are such that several thunder 
Cells combine to a single, giant cell. 

Typical Structure : Depends on terrain and meteoro- 

logical situation, 
life cycle: ^1 - 10 hr 

length; 50 - 150 km 
width; 'v 25 - 50 km 

speed of approach, severe gusts, and features are similar 
to those of multi-cell and large area thunderstorms. 

Special Features ; Large hail areas V7ith hail stones 
larger than 1.5” and tornadoes. Durations of 6 to 12 hrs 
(within multicell areas or isolated giant single cell) v75-th 
strong storm and hail periods from 1/4 to 3 hrs. 

7 9 4 5V 

2.1.4 , Thunderstorm (prefrontal) Squall Lines ‘ 

Occurrence : 80 to 500 km in advance of a cold front 

lifting particularly unstable warm and moist air mass or 
inversion. Cold outflov7s from developing line of thunder 
clouds reinforce cold air wedge far ahead of front by moving 
faster than upper level cloud portions . Resuiting line of 
rising columns of warm air masses in front of thunderstorm 
line and propagation of thunderstorm after thunderstorm, in 
sguall-line. 
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Typical Structure : Band of thunderstorm cells 

stretched perpendicular to direction of motion, 
life cycle; 1 r 10 hr 

length: 20 - 60 km in direction of motion 

v/idth: 250- 500 km (length of band) 

speed of approach: 15 - 60 km/hr . 

severe gusts: 20 - 30 m/sec at front of cold air 

outflow under and at leading edge of cloud line (roll cloud) 
Severe updrafts and dov7ndrafts similar to or more intense 
than those for thunderstorm cells (2.1.1). 

Strongest v;inds and V7ind shifts: similar to, or more intense 

than, those in and around 
storm cells of multi-cell 
thunderstorms (2.1.1 - 2.1.2), 
particularly in front, below, 

■ and just above the roll 

cloud near the lower leading 
edge of the squall-line clouds 

Sp ecial Features : Strong hail shov7ers V7ith strong 

down-drafts; tornadoes and funnel clouds are usual in strong 

line-squall thunderstorms . 

2.2 Tornadoes 

In order to distinguish operationally tornadoes and 
often less dangerous cold air funnel clouds and cold air 
v/ater spouts fron\ usually harmless dust devils, etc., v;e 
define them as strong vortices which interconnect from 
clouds to the ground or V7hich are connected to clouds with- 
out touchdown. Funnel clouds are included which form in 
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warm air sectors in contact with cold air masses. For 
cold air funnel clouds , see 2.3 belov/. 

Occurrence : Heavy multi-cell, area, and line squall 

thunderstorms (including supercell storms) , particularly, 
if rotating thunderstorm cells develop. Twisting updrafts 
and dovmdrafts next to each other in thunderclouds , parti- 
cularly at dry cold fronts and in regions of turbulent mix- 
ing of cold air with v;arm moist air at the cloud level in 
thunderstorm line-squalls. Require unstable temperature • 

profile with unstcible moisture prof ile' at lov7 altitudes. 

Typical Structure (See Figure 2.3.): 

path-length: 0.5 - 15 km (150 to 500 km observed on very 

rare occasions) . ■ 

path-width: 10 - 300 m (0,5 to 10 krn observed on very rare 

occasions v;ith multiple touchdowns) . 

speed .of app.roach: 15 - 30 km/hr (50 to 80 km/hr on rare 

occasions in fast-moving thunderstorms) . 

V7ind speeds: 20 - 70 m/sec (100 to 200 m/sec observed on - 

rare occasions) . 

updrafts: In suction areas and suction vortices. Velocities 

not really knovjn but may be estimated as 30 to 
100 m/sec (debris carried aloft can reach altitudes 
of 8 to 10 km in and above parent cloud) . Debris 
can be scattered as far as 10 to 20 km (50 km on 
rare occasions) by usually strong altitude v/inds 
of up to 200 km/hr at the top of parent clouds 
(8 to 10 km height). 

Speoial Features : Usually condensation and electricai 

discharge in center funnel and secondary suction vortices. 

Often , rotating thunderstorms can lead to cirea damage similar 
to that of tornadoes and cyclones by producing several tracks 



of suction marks. Tornadoes often occur in families in 
intervals of about 30 to 60 min; the individual families can 
exist for 20 to 50 min v/ith 3 to 8 tornadoes in a family. 
Over areas of water, tornadoes become tornadic v;ater-spouts . 

2,3, Funnel Clouds (Gold Air Funnel Clouds) 

Funnel clouds occur primarily in cold air masses 
(cold air funnel clouds) and v;arm moist air masses at line- 
squalls and cold fronts. The latter funnel clouds are 
typical precursors of tornadoes, and are therefore not 
included here (See 2.2 above,). Cold air funnel Clouds have 
no tornado potential and are usually low damage-producing - 
if damage occurs at all. They are: included here because of 
lack of study. {Laser velocimetry studies should be most 
helpful since rada.r signatures are usually absent) . 

Occurrence : Unstable moisture profiles below cold 

and dry air v;ith stable temperature lapse rate. Cold air 
circulation of cold core lov; pressure centers near weak 
upper-level troughs. Near shore lines if on-shore wind 
supplies moisture belov7 cold air masses. 

Typical Structure : Usually, cold air funnel clouds 

descend from parent clouds and rarely touch the ground. 
Diameters are typically 10 to 20 m at the base of the. parent 
cloud and taper to less than 3m near the ground. If touch- 
down occurs , there, is usually only very minor damage v;hich ’ 
is limited in area. Funnels do not nedessariiy extend 
vertically from the parent cloud. 


path-length: 1-4 km (up to 10 km observed on rare 

occasions) 

path-v;idth : 1 -5m (up to 20 m observed on rare 

occasions) 

speed of approach: 15 to 50 km/hr 

v;ind speeds: not really knov7n; may be estimated to 2 - 10 

m/sec because of minor damage. 

* Spacia'l Features : The duration of cold air funnel 

clouds is limited from 30 sec to 10 min (up to 15 - 20 rain 
observed on rare occasions) , Over areas vzith v;ater, cold 
air funnel clouds become cold air v?ater spouts. Cloudy 
associated v;ith funnels reach altitudes from 4000 to 8000 m. 
The funnels extend usually only 60 to 250 m below' the base 
.of the parent cloud (funnels of more than 1 km length have 
been observed on rare occasions). 

2 . 4 Dust Devils ^ ^ ^ 

Although most dust devils are of small scale and 
of small intensity, some can reach considerable strengths 
so that they can turn over cars and cause significant . 
structural damage. 

Occurrence : Very unstable temperature lapse rates 

near the ground in hot and dry air (strong near-ground 
inversion). 

Typical Structure : Hot-air, thermal vortices rotating 

clockvzise or counter-clockv/ise v/ith strong upv/ard convection; 
dust devils originate from the ground upv7cirds. 


path length: 0-5 km (up to 10 km on rare occasions) 

path v/idth; 10 - 100 m (up to 500 ra on rare occasions) 

height: 30 - 200 ra (very strong ones up to 1000 m v/ith 

invisible disturbances up to 200 m) . 

speed of approach:- 0 - 20 kra/hr (most dust devils move 

v;ith environmental v;ind if it is 
^ 5 km/hr); stationary at hill ridges. 

wind speeds: Not really knovm; may be estimated to 

5-30 m/sec in view of possible damage 

(up to 50 - 100 m/sec on very rare occasions). 

updrafts: Not really known; may be estimated at 10 - 50 

m/sec (up to 100 m/sec on rare occasions) in 
vi^W of possible damages . 

Special Features ; Most dust devils have life spans 
of 2 to 5 minutes; hov/ever, some may last for 30 min or more 
(up to 4 hrs on rare occasions). Dust devils do not show 
electrical discharges inside or at the base of the vortices. 
Usually, if discharge is observed, the phenomena is called 
a tornado. Often smaller updraft vortices seem to appear 
within the main vortex of dust devils. 

( 818 ) 

2.5. Vortices Generated by Density Currents 

Strong vortices, like tornadoes, are often associated 
with gravity or density flovrs in the atmosphere which are- 
generated by severe temperature differences. Usually, the 
necessary flov;s result from the cold outflov/s of lax"ge area 
thunderstorms and line squalls and from cold air (katabatic 
winds , Chinooks , Santa Anas , boras) , particularly v;hen 
supported by steep dov/n slopes in , the terrain* Similar 
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types of strong vortices are produced, by density flows, 
for example, in airborne avalanches of pov7der-snow, avalanches 
of liquid and solid volcanic materials that evolve gases 
continuously (nuee ardente) , and large sand and dust . 
desert storms (haboob) . Since these phenomena produce, 
high velocity, vortex-like circulations similar to those 
observed in severe thunderstorms , they are briefly outlined 
here and their common features are described below. 

(8 X3 ) 

2 . 5 . 1 . Airborne Powder-Snow Avalanches ^ * ■ 

Dense suspension of air surrounded snov; particles 
moves like a heavy gas with very little internal friction. 
Surrounding warmer air and irregularities of terrain cause' 
friction along outer surfaces. The result is considerable 
turbulence and increase in the speed of the central core 
of the cold, high-density flow, v;ith spin-off of both, 
vertical and horizontal vortices. . . 

bulk-flow velocity: 30 - 100 m/sec 

vortex-flov; velocity: 2 x bulk velocity (up to 200 m/sec) 

vortex-flow duration: 2 •- 3 sec for each cell 

vortex cell size: 0.5 - 30 m (estimate) 

2.5.2 Nuee ardente ''^ (7,11) 

Avalanche of heavy volcanic material ejected as a 
viscous-, - auto-explosive liquid that continuously evolves 

* I ~ ^ ^ 

*''Nuee a.rdcnte" from french : nuee - “storm'* , and cxrdent - 

"fiery, red," “burning"; thus "firey storms" (not to be con- 
fused V7ith "fire storn'*" : strong surface st or tas genex'atcd 

bv,,,the_.suction of l^irge area fires) . ' 

— ■' : ; { ■ 'V ' ■ I - I ", r- 


gases . Emitted gases practically eliminate contact friction 

of the particles and cause avalanche to spread almost as 

freely as a true liquid. Irregularities in density and 

terrain result in considerable turbulence and spin-off of 

horizontal vortices. 

bulk-flow velocity: 3-10 m/sec 

vortex-flov; velocity: 70 - 100 m/sec 

vqrtex-flow duration: not really known (-1 - 20 sec, 

estimate) 

vortex cell size: not- really known (0.5 - 30 m, estimate) 

front of avalanche cloud: many short 0.‘5 - 30 m ?) 

horizontally-rotating vortex 
. cylinders v:ith axes of no two 
cylinders parallel. - 

2.5.3* ' Sand ‘and Dust Storms (haboob) ^ . 

Huge sand and dust storms are usually generated by 
cold air outf lo'ws of decaying thunderstorms during desert 
seasons of high convective instability, or by cold air, kata 
batic winds (Chinooks, Santa Anas, or boras) resulting 
from large high pressure areas in northern or so\athern 
cold regions of the globe. They may stretch over v/idths 
of more than 150 km and can rise to the level of the 
cloud base (condensation base) . A. in the case of nuees 
cirdentes and airborne powder snow f.v' ..anches, a dense 
suspension of air-sircrounded sand and dust pjirticles evolves 
vrhich moves like a heavy gas vjith very littlo internal 

''Haboob" froje- Arabia i hebbe - "blovm". ' 
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friction. Surrounding various air masses and air pockets 
as v;eil as irregulcirities in the terrain cause friction 
areas v/ith resulting build-up of intense vortices, parti- 
cularly at the turbulent fronts of these systems. Often, 
"eddy tornadoes" form in the lee of local minor mountain 
system.s. 

bulk flow velocity: 5-20 m/sec (or larger, depending on 

the velocity of the cold air outflov?) . 

vortex-flow velocity: 20-100 m/sec (or larger on rare 

occasions) . 

vortex-flow duration: not really known (~ 5 to 60 min. for 

tornado-like structures estimated) . 

vortex cell size: - 10 - 300 m (or larger on rare occasions 

for tornado-like structures) . 

2.5.4 Typical Vortex Structures in Density Currents ^ ^ ^ 
Atmospheric vortices caused by density currents 
resulting from temperature differences are similar to those 
observed in true liquid systems at the interfaces of gravity 
flows, in the following, three structures are briefly 
described V7hich are typical for such flows. They are intend 
ed to provide with a general and qualitative idea about hov; 
and v/here intense vortices can be formed in the atmosphere.. 

'i C A T r. Tl T (8,14,15) 

2.5. 4.1. Surge Heads 

Heavy cold air flow into X'/armer air xnasses can pro- 
duce surge heads yzhich develop characteristic protruding 
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Figure 2.4a. Side View of Cold Air Surge Head Moving 
from Left to Right. . 



Figure 2.4b. Top View of Cold Air Surge Head Moving 
from Left to Right. 
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lobes separated by incisive undercutting clefts and warm 
air tunnels (see Figs. 2.4a and 2.4b) along the ground 
surface. Eddies with vertical axes often form v/ithin the 
clefts. The surge head may become a system of oppositely 
rotating longitudinal vortices which are practically 
identical to the internal structure of the nuee ardente. 
Similar atmospheric vortices can also be caused by sea 
breeze fronts, by dust walls in dust storms, and by the 
leading edges of steep cold fronts. The resulting phenomena 
may be strong dust devils, funnel clouds, and tornado- 
like Structures, either occurring isolated pr in rov?s along 
the front of the cold air masses. Tornado-like vortices 
can develop in the presence of as little as one percent 
density difference betv?een the cold and v;arm air masses, 
corresponding to temperature differences of only 3°C, 
Temperature differences of more than 13°C (corresponding to 
density change of > 4%) have been observed at the front of 
large haboobs fed by cold air outflovjs from large thunder- 
storms. 

2.5 . 4 . 2. Ri s ing Tunne Is ^ ^ ^ ^ ^ 

Vertical vortices in density currents can occur when 
the warm air entrained in the tunnels of surge head rises 
(see Fig. 2.5). This rise can be caused by eddies in the 
wake of an obstruction- in the path of the cold air outflow 
and can result in giant dust-devil type vortices t If the 
volume of engulfed v/arm ciir is large enough , the vortex 
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Figure 2. 5 . 


Plow (from left to right) of a Surge Head 'Over 
a Topographic Obstruction. (Eddies in the wake 
of the obstruction can cause warm air tunnels 
to rise like a giant dust devil.) 
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Figure 2.6 


Vortex Pair Formed in a Density Current by Inter- 
action vaith a Topogrciphic Disturbance. (Vortex 
tubes may originate in the upper or lower 
circulation of the surge head .shov/n in Ficfure 2.4a.) 






can penetrate the top of the surge head. Then v/arm air 
pushed up by the leading edge of the cold air outflo;*; may 
be added to the vortex so that it can reach the condensation, 
base. Here, the vortex can gain further energy by con- 
densation of moisture to nev; cloud elements. At this stage, 
the vortex has acquired the characteristics of a full- 
fledged tornado. Although it can occur everywhere, this 
mechanism is probably more likely to operate in the presence 
of topographic disturbances. Tunnels may rise by buoyancy. 
The turbulence inherent in the surge head can impart easily 
a circulation pattern to the acending convection current 
and, once stretched, the vortices can intensify considerably, 

2. 5. 4. 3. Vortex Pairs 

Vertically rising vortex pairs can be formed v;hen 
horizontal vortex tubes are lifted out of the local boundary 
layer by a topographic disturbance (see Fig. 2.6), The 
horizontal vortex tubes may originate either in the upper 
or in the lower lobe-circuilation of the surge head associated 
with, a cold air outflovz as shown in Fig, 2.4a, 

2. 5. 4. 4. Vortices from Clouds ^’^^^^ 

Vertical vortices reaching dovm from clouds can be 
formed by similar mechanisms as those described above if 
strong rotating dov/ndr aft s (vertical doxmward surge heads?) 
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engulf warm updrafts with the same axis and sense of 

(9) 

rotation (see Fig. 2.7). This situation may appear in 
rotating thunderstorms and vrhen strong vertical cold-air 
outflows from huge thunderstorms and line-squall thunder- 
storms are accelerated by hail, snow and rain. Usually, 
tornadoes are the result at the base of these regions 
where cold downdrafts and v/arm updrafts mix in a swirling 
motion. 

NOTE : Valuable criteria have been quoted in ref. 20 

for given buoyant forces with access to’’ large scale vortici- 
ties to develop into concentrated vorticities; namely, for 

.15: buoyancy dominant - no significant vortex 

generation, 

>0.9: no concentration of vorticity by vertical 

motions - vzhole flo'W sv/irls, 

0,15 < U^/U_ 0.9: possible generation of concentrated 

vorticity. 

I ■ 

Here/ = tangential component of horizontal velocity of 
environmental air and - characteristic vertical speed of 
buoyant air mass . . 

Considerable more observatiGns are required of the 
vortices associated with frontal boundaries of density flov?s 
if a reasonably confident understanding of their true nature 
and dynamics is to be achieved. Typical guestions arer Do 
they occur in front of or. behind the surge heads? Is their 
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Figure 2. 7.- Possible Mechanism for Tornado Formation by 
Co-^Axially Rotating Warm Updrafts and Cold 
Downdrafts* 


A-II-31 


rotation cyclonic or anticydonic? Do they have dynamic 
connections v;ith the cloud base and/or the air near the 
ground? Are they favored downv/ind of significant topo- 
graphic disturbances? Are twin-funnels very common? Are 
smaller subsidiary vortices often associated vrith large 
dust devils? . 



3. TYPICAL STRONG WIND SHIFTS AND GUSTS ALOFT AND NEAR 
THE GROUND 

Large local changes in the direction and the speed 
of^v/inds aloft and near the ground occur primarily by two 
mechanisms: (a) high velocity, veered (turned clockwise) 

v/inds brought down by prevailing vortex loops or vortices 
in the lee of topographic disturbances and (b) vortices 
and turbulence generated by v/ind shears and by topographic 
obstructions in the main stream of strong gradient, v^arm 
katabatic, jet stream, and hurricane winds. All other dis- 
turbances associated v;ith the mixing phenomena occurring 
between warm and cold air masses and pockets fall under 
■disturbances in density flov/s (see above). Gusts are related 
to turbulence in air flov/s, and v;e consider here moderate 
to extreme turbulence v/hich is usually associated with v?ind 
speeds of 10 m/sec (20 kts) and larger. 

Typical gust factors are shown in Table 2 as functions 
of wind speed and average time periods used for obtaining 
data. From Table 2, it is apparent that, for average time 
periods between 1/2 and 30 sec, the typical gust velocities 
are given by 



gust 


i (1.44 ±0.11) 


These values app3.y for the horizontal wind component in 
gusts. ■ 
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Time Period; sec 



1.50 


1.47 


1.48 


1.61 


1.49 


1.4 8 


1.48 



1.35 


1.37 


1.35 


1.36 


1.35 


1.42 


1.36 


1.36 



TABLE 2. Gust Factors Averaged Over Short Time Periods 
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The vertical velocity component follows the same relation 

for altitudes - 30 m. Decreases of the gust factor v;ith 

increasing altitude are not considered since they apply 

only for ground friction eddies in stable air streams. 

High speed horizontal gusts have typical durations of 1/2 

to 30 sec; the corresponding vertical durations are 

- 1/3 -times the horizontal periods . Typical distances 

over which the average v/ind speeds increase to peak gust 

velocities are between 15 and 70 m and correspond to v;ind 

shears betv;een 1/10 and 3 sec Table 2 can be used to 

estimate typical peak velocities in gusts for gi'^en wind 

speeds. Thus, for Santa Ana v;inds of 30 - 35 m/sec,* down- 

draft gusts invalleys ' of up to 40 m/sec are possible and 

have been observed. Typical directional changes of v/inds 

( 2 ) 

in gusts are 15® to 30® near the ground . But changes 
of more than 160° have been observed under roll clouds 

(2) 

caused by the rotor flow in the V7ake of dense-bamer flov7 . 

. (2) ■ 

3.1. Typical XVave Streamxng and Rouor Streaming 

Wave streaming and rotor streaming are vortex- 
shedding phenomena and occur usually in the lee of large 
topographic disturbances (see Figs. 2.1a and 3.1b). In 
the case of rotor streaming , counter-flov; develops at 
the lee side of the obstruction near the ground, and 
the f lovi is more turbulent than that below the rotor 
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Altocumulus Lenticularis 



Figure 3.1a. Typical Wave Streaming' in the Lee of a Large 
Topographic Disturbance. 


Rotor Cloud 



Figure 3i’lb. Typiccil Rotor Streaming in the I.ee of a Large 
Topographical Disturbance (vortex shedding) . 
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clouds in wavestreaming. 


wind shifts: from 25® to 160®, particularly under rotor 

clouds , 

velocities in gusts : depend on the average v/ind velocities 

(See Table 2.)* Note: momentary 

calms can occur below rotor clouds and 
roll clouds, as v/ell as in the near 
lee of mountains v;ith subsequent wind 
shifts of -90® to 135® to the direction 
of the previous gust. 

Typical Features : Severe downdrafts can occur momen- 

tarily and periodically below and near the front of roll- 
clouds V7ith velocities on the order of (and up to 2~times) 
that at the mountain ridge. Directional shifts near and 
at the ground are usually large in the near lee of the 
obstruction and under the rotor and roll clouds. P7ind 
shifts of 90® to 135° have been observed at the ground 
between observation stations V7hich v;ere 800 m apart; for 
example, at an airport in the lee of a mountain range.' ' 

3,2,^ Turbulence in the Lee of Mountains 

Moderate Turbulence : Found for winds of < 12 to 25 m/sec at 

the ridge of mountains at distances of up to 250 to 500 km 
leev/ard, from the ground upi to the tropopause and- as much 
as 3000 m above the tropopause. 

Severe Turbulence ; Found for winds of >25 m/sec at the 
ridge of mountains at distances of up to 80 to 250 km leeiifard 
and from the ground up to the tropopause. 
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BULK MOTION 


100 - 2500 m. 


•^1 /J 


Typical Circulation in Medium Gust Cells 
as reviev;ed in a coordinate system moving 
v;ith the circulation cells 



100“2500 m. 



GUST TOE 


Figure 3.2b. 


Typical Section through Medium Gust. Colli 
(horizontal motions near Coll Bou: laries) 




Extreme Turbulence ; Pound for winds of > 25 m/sec at the 
ridge of mountains at lov/ levels in the lee and in and near 
rotor clouds. Sometimes severe turbulence is found v?ithin 
a layer betvreen the height of the jet-stream core/ • 2000 'm' : 
below the core and about SO to 250 km tov;ards the cyclonic 
(cold) side of the jet core. 

3.3. Circulation Gust Cells ^ ' 

Circulation gust cells are caused by frictional flov; 
near and at the ground. Usually/ higher sjjeed wind is 
brought to the surface more or less periodically from aloft 
(see Figs. 3.2a and 3.2b). Usually high-speed gusts are 
brought dovm behind the leading edge of the gust fronts. 

wind shifts; In the northern hemisphere/ clockv;ise (veered) 
V7ith increasing speeds, counter-clockv/ise 
(backed) v;ith decreasing wind speeds; 
typically 15° to 30°/ but up to 45° observed 
on some occasions. ■ 

v;ind speeds : See Table 2 for gust velocities as a function 

of vzind velocity. 

average size; 100 - 2000 m for small to medium size eddies 
. of strong gusts . 

Special Features : For medium size gust cells, the 

v/hole v/ind field, is rather variable in both speed and dir- 
ection. If the upvzard convection in gust cells extends to 
large altitudes/ major gusts with tv/ice the average v/ind 
spaed can occur near and at the ground in large eddies over 
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~ 30 min average recurrence intervals. The typical life 
times of gust cells are not v;ell knov/n, but may be from 
5 to 50 sec. The turbulence usually associated with 
gusts can extend to altitudes of 1000 ~ 2000 m. 
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20. Abstract 


scanned for velocity (and possibly turbulence ) , and is intended 
for diagnostics of the lower altitude region below the storm 
and in the out-flow region. 

The essential conclusion of the study is that aircraft flight 
velocities are high enough and severe storm lifetimes are long 
enough that it is possible for a single airborne Doppler system, 
operating at a range of less than about 20 km, to view the storm 
area from two or more substantially different aspects before 
the storm characteristics change appreciably. This permits 
resolution of the spatial distribution of the velocity vector 
by a single coaxial system operating in a fly-by mode. 

In addition, it is concluded that when either or both peak 
power rather than mean pov/er and/or peak data handling rate 
rather than mean rate are limiting 'factors, the optimum signal- 
to-noise as well as optimum velocity resolution are to be obtained 
by using pulse lengths comparable to the overall dimensions 
of the viewed region and, as v/ell, by operating in a continuous 
azimuthal scan mode* Range resolution is to be recovered in 
this mode by post-flight image reconstruction ana.lysis. 
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